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Carbon-14 Carboxy-Labeled Polysaccharides’ 
Joseph D. Moyer and Horace S. Isbell 


Procedures are given for the preparation of C' carboxydextran and C™ carboxyvinulin. 
These materials have C-labeled carboxyl groups in place of the reducing end groups 


present in the parent polysaccharides. 


The substances can be prepared cheaply and should 


find many applications as tracers in biological and chemical research. 


1. Introduction 


A process for labeling polysaccharides with carbon- 
14 was described briefly in a prior note [1] * and in a 
patent assigned to the | nited States Secretary of 
Commerce [2]. Subsequently, the labeled products 
have been furnished to research workers in other 
laboratories, but the procedures for making the com- 
younds have not been deseribed in detail. 

The present paper reports the methods used for 
the production of C'™ carboxvdextran and of C™ 
earboxvinulin. These materials have already found 
yse in the study of adsorption on blood platelets and 
cells, as well as for the study of kidney function [3, 4]. 
Other applications can be envisioned and will be 
realized in the future. 

In the labeling process, reducing end groups pres- 
ent in the poly saccharide combine with C-labeled 
évanide; the resulting cvanohydrin, on hydrolysis, 
yields a polysaccharide having a radioactive carboxyl 
group. The activity of the labeled product depends 
on the equivalent weight of the polysaccharide and 
on the activity of the cvanide used in the synthesis. 
The radioactive product differs from the parent poly- 
saccharide in that it has a C-labeled carboxyl] 
goup in place of the reducing end group. The 
presence of a single carboxyl group in a large mole- 
cule has little effect on the biological behavior of the 
substance. Hence the carboxy-labeled polvsaccha- 
rides can be used as tracers in biological and chemical 
research. Prior to starting the svnthesis, the com- 
bining weight of the polysaccharide is determined by 
reaction with radioactive cvanide by the method 
already described [5]. The polysaccharide is then 
treated with C' labeled cvanide and, after removal 
of the excess cvanide, the product is isolated. The 
methods used for purification and isolation of the 
product depend on the properties of the material 
being labeled. The procedures described here can 
be used, with small changes, for almost any water- 
soluble polysaccharide. 


2. Experimental Procedures 
2.1. Materials 


Carbon-14-labeled evanide was prepared by the 
previously described method 16). The solution used 
contained 0.065 mM _ of sodium evanide plus 0.15 
mM of sodium hydroxide per milliliter and had an 
activity of 1 me/ml. 
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The ion-exchange resins used for release of hydro- 
gen cyanide from the reaction mixtures were condi- 
tioned by successive treatment with 5-percent aque- 
ous sodium hydroxide, 5-percent aqueous sodium 
evanide, water, 5-percent aqueous hydrochloric acid, 
and water until the wash liquor was neutral and free 
from chloride. 

The clinical dextran was supplied by the Commit- 
tee on Shock, National Research Council. The sam- 
ple, designated NRC-3, was part of a supply made 
for experimental studies by the Commercial Solvents 
Corporation. 

The inulin used was a commercial product. 
labeling, it was recrystallized from water [7]. 


Before 


2.2. Determination of Cyanide-Combining Power 


Approximately 5 mg of the lyophilized polysac- 
charide was placed in a weighed test tube which had 
been constricted near the middle for sealing in a 
flame. The sample and tube were dried at 80°C 
under a pressure of less than 0.1 mm, and the weight 
of the sample was ascertained. The analyses were 
set up in triplicate, with controls containing 0.02 
mg of p-glucose and with blank determinations hav- 
ing like quantities of the reagents only. 

A buffered sodium eyanide-C™ solution was pre- 
pared from the 0.005 N carbon-14-labeled cyanide 
solution by adding crystalline sodium bicarbonate 
equivalent to the excess sodium hydroxide. A three- 
to five-fold excess of the evanide solution (approxi- 
mately 0.1 ml of the solution, containing about 0.2 
uc of carbon—14) and 1 drop of toluene were added 
to each tube. The tube was sealed and stored at 
room temperature. After 10 days the tube was 
opened, a few drops of 10 percent aqueous formic 
acid was added, and the solution was concentrated 
in a current of air. The residue was dissolved in a 
few drops of aqueous formic acid, and the solution 
was frozen and lyophilized to remove the last possible 
trace of hydrogen cyanide. The material in the tube 
was dissolved in formamide, and the radioactivity 
was measured with an internal, gas-flow, propor- 
tional counter [8]. The equivalent weight was ob- 
tained from the relationship: 


. : ‘ (cps) VW 
Equivalent weight ——. 
; (CPs) w-G 


where (CDS)e@ is the net count of the control contain- 
ing G millimoles of p-glucose and (eps) is the net 
count for the polysaccharide sample of weight W. 








2.3. Carbon-14 Carboxy-Labeled Carboxydextran 


Ten milliliters of a 0.065 M sodium cyanide-C"™ 
solution was placed in a 100-ml flask containing 2.6 g 
of lyophilized dextran (NRC-3). The flask was 
sealed by a glass stopper and a ring of silicone grease. 
The dextran was dissolved and stored at room tem- 
perature. 

After 4 days, a 10-yl sample was withdrawn for 
analysis, and the flask was immediately restoppered. 
The sample was transferred, together with washings 
of the pipet, to a small test tube and acidified with 
3 drops of 10-percent formic acid. ‘The mixture was 
evaporated to dry ness in a jet of air. ‘To insure com- 
plete removal of hydrogen cyanide, the residue was 
twice dissolved in 0.5 ml of water and reevaporated 
to dryness. Radioassay of the residue, dissolved in 
1 ml of formamide, showed that approximately 
0.67 me of carbon—14 had been fixed. This was 88 
percent of the amount expected (0.76 me) from the 
specific activity of the evanide (15.2 me/mM) and 
from the cyanide-combining power determined by 
the method of section 2.2. 

The reaction mixture in the 100-ml flask 
allowed to stand an additional 3 days to assure com- 
pletion of the reaction. The liquid was then frozen 
by immersion of the flask in a liquid-nitrogen bath 
with gentle swirling. The stopper was removed, 
5 ml of conditioned Amberlite [R-120(H was 
added and the flask, A, was connected to a receiving 
flask, B, as shown in figure 1. The receiver B con- 
tained excess sodium hvdroxide solution (1 mil-of .V 
NaOH). Flasks A and B were immersed in liquid 
nitrogen and the system was evacuated with an effi- 
cient oil pump. Stopeocks E and C were closed, the 
liquid-nitrogen bath was removed from A, and the 
contents of the flask were thawed. After thorough 
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mixing, the contents were again frozen in liquid Nitro 
gen. Stopeock E was then opened and the liquid 
nitrogen bath was removed from A, but not from R 
As the frozen mixture became warmer, water gyj 
hvdrogen evanide were transferred to B by lvophil. 
zation. After the lyophilization was complete (5 hy)| 
stopcock C was closed, and 10 ml of water was added 
to A (through funnel D) and mixed with the conteyy 
of the flask. The liquid was frozen in liquid Nitrogen 
the system was evacuated, and the Lyophilizatigy 
process was repeated. Finally, flask A was remoya 
from the apparatus, and flask B, containing the! 
recovered cyanide, was stoppered and set aside for 
future use. 

The contents of A were passed through a colump 
containing 10 ml of a 1:1 mixture of Amberlit. 
IR-120(H) resin and Duolite A-4 resin.’ The flask 


and resin were washed with water. The Solution was 
acidified with a drop of 10-percent acetic acid ani 
lvophilized. The residue was dissolved in 50 ml 
water and again lyophilized The second lvophili.’ 
zation was to remove the last trace of hv droge 

evanide.) The product weighed 2.29 ¢ and had, 

activity of 0.28 we/meg. 


2.4. Carbon-14 Carboxy-Labeled Carboxyinuliy | 


Seven grams of finely powdered, lyophilized, y 
crystallized inulin was added slowly with constal 
swirling to 43 ml of an ice-cold cyanide solution in 
100-ml flask. The evanide solution contained 04 
mM of sodium cyanide—C™ (7.2 me) and 0.3 mM ¢ 
sodium hydroxide. The flask was capped witha 
adapter carrying a stopcock for evacuation. Thy! 
contents were frozen in a liquid-nitrogen bath a 
the flask was evacuated 
and the flask was allowed to come to room temper 
ture. It was then heated in a water bath, wit 
constant swirling, until the inulin had dissolve 
After two days at room temperature, the inulin th 
had crystallized from solution was redissolved bj 
warming. This was repeated every other day f 
Sdavs. At this time, an analysis on a 20-al samp 
showed that 5.2 me of carbon—14 had been fixe 
This amount corresponds to approximately the activ 
itv anticipated from the evanide-combining power: 
inulin (previously measured) and the specific activi 
of the evanide used, 

The contents of the flask were frozen in liqui 
nitrogen and, after the addition of 5 ml of cond 


The stopcock was close 


tioned Amberlite IRC 50 H resiti, the flask wir 
connected to a second flask (B of fig. 1 containing 
| ml of N NaOH. Flasks A and B were immers 


in liquid nitrogen, the svstem was evacuated, and! 
hvdrogen evanide liberated in A was transferred t 
Bb by the procedure described in connection wit! tl 
preparation of labeled dextran. After removal of a 
evanide, the residue in A was extracted at 40° tod 
(' with portions of water totaling 150 ml. Ti 
solution was filtered and passed through a colun 
containing 5 ml of a 1:1 mixture of Amberlite IR 
50(H) resin and Duolite A—4 resin. The combin 
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quid Nitro. fluent and washing were lyophilized after the addi- | 
: The residue was | 


the liquid: “ion of 3 drops of acetic acid, 
Ot from R dissolved in 100 ml of warm water. A 20-yl aliquot 


Water and) was removed and counted for radioactivity in l ml of 


'v lyophili. formamide; the remainder of the solution was 
ete (5 hr) | jyophilized. The residue weighed 7.1 ¢ and _ its 
> ~ 


Was added radioactivity showed the presence of 5.14 mc. 

1€ CONtents| Elution of the resin column with 5 ml of 5 percent 
d nitrogen HCl, and washing, resulted in the recovery of 0.16 
Philization me of activity. The cyanide collected in flask B had 
iS Temoved an activits of 1.08 me. 
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thank Nathan L. Drake of the University of Mary- 
land for helpful advice. 
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Heat of Formation of Sodium Calcium Aluminate 


Edwin S. Newman 


The compounds 3CaQ-Al,O0O3, NasSO, V 


(Thenardite), 


Na.O-8Ca0-3Al,03, and CaSO,, 


snhydrite, were prepared, and their heats of solution in HCl, 26.61H,O were determined. 


The heat of solution of Na.SO, V 


in water to form Na.SQ,, 1000H,O was also measured. 


From these and other data in the literature the heats of formation of Na,O-8Ca0O-3Al,0; 


were calculated to be 


and NaSO, V 
1. Introduction 


Investigations of the alkali phases in portland ce- 
ment clinker, summarized by Newkirk [1],' have 
hown that the compound Na,O-8Ca0-3Al,O,, here- 
after abbreviated NC,A,,’ is formed. As part of a 
ntinuing investigation of the thermochemical prop- 
erties of compounds occurring in hydraulic cements 
and their hydration products, the heat of formation 
of this compound was determined. No prior ther- 
mochemical data on NC A, have been found. Meas- 
ements were made at 25° C of AH of the reaction 
3C,A(e) 4 
NasQO,(c. V). 


VC.A,(c) + CasO,(c, anhydrite) 


(1) 


This reaction occurs spontaneously [4] at high tem- 
peratures, where, however, the high-temperature 
form of calcium sulfate is present and liquid sodium 
sulfate is formed. The value of AH, was obtained 
wv measuring the heats of solution of the four sub- 
dances in hydrochloric acid solution. The heat of 
formation of NC A, was calculated from this value 
of AH, and the heats of formation of the other three 
substances taken from the literature. 


2. Materials, Apparatus, and Procedure 


Tricalcium aluminate was prepared by the method 
sed by Thorvaldson, Brown, and Peaker [2]. The 
gnition loss of this material after grinding was 0.08 
percent, assumed to be water in the compound 
(AH, formed by reaction of the anhydrous alumi- 
ate with the moisture in the air. Thenardite, 
\aSO, V [3], was crystallized from a 50-percent so- 
ition of Na,SO, at 50° C by the slow addition of 
ot alcohol. The mixture was left on the steam 
ath for 2 hr and then filtered. The crystals were 
lned overnight at 100° C and eround to pass a No. 
They were essentially free from inclu- 
ons, as judged by microscopical examination, and 
the loss in weight on heating for 2's hr at 950° C was 
(2 percent. In the measurements and calculations, 
is material was considered to pure Na SQ,. 
‘he compound NCA, was prepared as described by 
Newkirk [4], except that after heating once at 1,200° 
Uthe preparation was autoclaved for several hours 
i 300 psi before repeated heating at 1,300° C with 
itervening grinding and sieving. Microscopical 
tamination and chemical analysis indicated that 
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2576 and 


| surface of a water bath controlled at 
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330.92 keal/mole, respectively. 


the product was homogenous and of the proper 
composition. Anhydrite was prepared by first 
precipitating gypsum from a hot solution of CaSO, 
prepared from the hemihydrate and then converting 
the gypsum to anhydrite by heating for 2 hr at 
1,000° © [5). 

Details of the calorimetric apparatus and proce- 
dure have been given elsewhere [6]. The calorimeter 
used in this work was a 1-pt wide-mouthed vacuum 
flask cemented to a metal flange. A brass cover 
with tubes for stirrer, thermometer, heater, and in- 
troduction of the sample was bolted to this flange, 
and the calorimeter was submerged 4 in. below the 

25.0° C. The 
cover was made watertight by means of an O-ring 
seal. The temperature of the bath remained con- 
stant to better than +0.002° C during a determi- 
nation. The calorimeter charge was 425.0 g of 
HCl, 26.61H,O (2.000 normal at 25° C). The start- 
ing conditions were so adjusted that the final tem- 
perature of the calorimeter was as close as possible 
to 25.0° C. This procedure was adopted both to 
avoid the need for correcting the measured heat of 
solution to the standard temperature of 25° C and 
to prevent the distillation of water from a warmer 
surface of the acid solution to condense on a cooler 
inner surface of the metal cover. Because there 
was no other free water surface within the flask, 
condensation in the reverse direction was unlikely. 

The weights of the samples for the heat-of-solution 
determinations were in the proportions required by 
the stoichiometry of the equation for a’weight of 1.42 
g of C,A. First the NC.A; was dissolved in the calo- 
rimeter and the final rate of temperature rise deter- 
mined to complete the measurement of its heat of 
solution. This final rate became the initial rate for 
the CaSO, experiment. The sample of anhydrite was 
introduced and the determination was completed. 
This procedure was followed also with a new acid 
charge for the products on the right side of the 
equation. Thus the final calorimeter solution for the 
measurements with the reactants was of the same 
composition as that for the products, and no dilution 
experiments were necessary. The experimentally 
determined energy equivalent of the calorimeter and 
acid charge was used in computing the heats of solu- 
tion of NCA, and C,A. For CaSO, and Na.SQO, it 
was increased by 0.28 cal/deg, the estimated heat 
capacity of the solid previously added. Since this 
increase was less than 0.1 percent of the total energy 
equivalent, the error introduced by this approxima- 
tion was negligible. 











of HC1-26.61H.O 


TABLI ! Hleats of solution n AA0 @ | P 
NagO-SCaO0-3AbhOs(S16.51 CasO,4(136.15 $CaO-AleOs (270.20) Naps Og V (142,06) Th 
Sample number ~ | 
Sample Heat Sample Heat of Sample Heat Sample Heat of | 
wt t olution weight solution weight olutior weight solutior 
, } dia t 
| r 1 qd li; ¢ iid cal , 
i ] 1 4220) $3 0. 241 13.00 1. 4200 736. 43 0. 251 32.43 
i 2 1. 4233 733. 81 3S 12.88 1. 4189 735. 49 2518 33 0 
i 3 1. 4254 734.0 2385 12. 3¢ 1. 4220 735. 27 ti 39 
t 1. 4262 4.11 2378 12. 83 1. 421 735. W2 At 32. 6 
§ 1. 4276 4.2 2407 12. $2 1. 4228 735. V1 2544 32 35 { 
AV $3.0 12. 69 i S0 9 6 
S Dp. 0.1 0.14 0.19 P 
al cal mole il mole nole 
H 40 2S l m= +41) OLD 199. 11 0.051 1 (40 +0. 017 
Numbers in paren molect we or C.AH 


3. Heats of Solution 


The results of the measurements of the heats of 
solution in hydrochloric acid are given in table 1. 
The negative sign before values in this table indicates 


that the heat of solution was endothermic. The 
value shown there for C.A when corrected for the 
calculated C,AH, impurity becomes 736.9 cal/g. 


This value, which was used inh computation, ha be 
compared with some previously obtained results. 
In 1930, Thorvaldson, Brown, and Peaker [2| obtained 
736.0 eal/g® for the heat of solution of 2.86 ¢ of C.A 
in 600.4 ¢ of HCI, 20H.O at 20°C.) In 1939 Cirilli {7 
obtained 735.1 eal/ge for the heat of solution of C.A 
in N HCI at 25° CC. Recently Coughlin [8 
obtained 738.0 cal/g¢ for the heat of solution of 1.43 
gf of C.A in 1936.2 ¢ of 4.360 M HCLat 30° C. These 
four values, including the one obtained in this work, 
are in reasonably cood agreement 

No data were found in the literature on the heat of 
solution of sodium sulfate in’ hvdrochlorie 
Two determinations in this laboratory with larger 
samples (1.42 g) than were used for the data in table 
| gave values of 32.38 and SLS7 cal/g, 
tively. 

Newman and Wells obtained 12.75 eal/g¢ for 
the heat of solution of CaSO, in HCl, 26.61H.O, and 
Southard [9] obtained 12.28 eal/g in 2.08 N NCL. 
These values are in substantial agreement with those 
in table 1. 


Oo. 


acid. 


respec- 


4. Heats of Formation 


Krom the heats of solution in table 1, the heat for 
reaction | at 25° C is calculated to be AH SO 

0.19 keal/mole. To determine the heat of forma- 
tion of NC A, the heats of formation of the other 
three members of eq 1 must be known. 

The heat of formation of C.A civen S61 
keal/mole in Circular 500 of the National Bureau of 
Standards [10!. This value was in computa- 
tions. Although Coughlin [S8!} gave the value —850.04 

1.2: 853.46 1.23 keal/mole depending on 
the value used for the heat of combustion of calcium, 
3.18 keal, twice his 

of C.A from the 


Is is 


used 


> or 
these values appear to be low bry 
value for the heat of formation 
oxides. 
, 3 The 


value given in 2 by 1 ; slated 


her " 1.184 


76 


The heat of formation of NasSO, V_ was obtained 
by measuring its heat of solution in water to fory 
Na.SO,, LOOOH.LO, the heat of formation of whie i 

331.215 keal {10}. Five determinations gay 
295.0 2 9 cal/mole for the enthalpy Of solutioy 
25 Coughlin’s value at 30° © [11), —G289 

was calculated to 302.2 cal/mole g 
The heat of formation of NasSO, V was 
15+. 0.293 330.922 keal/mole 


at Ss; 
cal mole, 
25 +} 


be 
taken as —331.215 
The heat of formation of CaSO,, anhydrite, j 


Ty 


$42.42 keal/mole [10]. From these values for the acry 
heats of formation of CyA, Na,SO, V, and anhydrite; desi 
and the heat of reaction AH,, the heat of formatigy| NB] 
of NCA, from the elements at 25° © is calculated ty port 
be —2576 keal/mole. | resu 
stab 

5. Summary j whe 

sim 

The compounds 3( ‘aO-ALO,, NaSO, V_ (thenary of | 


dite Na ()-N¢ ‘af .SALO . and (‘ast d.. anhyvdrit D 
were prepared, and their heats of solution in Hj tion 
26.61 HO were determined The heat of solution of met! 
NaoSO, V in water to form Na.SO,,1000H,O was sam 
also measured. From these and other data in the first 
literature the heats of formation of Na,O-8CaQ) 30) 
3ALO, and Na.SO, V were calculated te be —257h weig 
and 330.92 keal/mole, respectively. The 
mas 
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«« —\Thermal Degradation of Polyacrylonitrile, Polybutadiene, 
and Copolymers of Butadiene With Acrylonitrile 
and Styrene 


Heat of 


Solution 


cal/g 
32. 43 
33. 02 

~32. 82 
32. 


82. 35 


S Obtained! 
er to form 
vf which is 
ions vay 
f solutioy 
L}, —6289 
al/mole 
( ), V Was! 
1 mole 

ivdrite, js! 
les for the} 
anh yvarit 
formatioy 
culated t 


Sidney Straus and Samuel L. Madorsky 


Polybutadiene and a copolymer consisting of 76.5 percent butadiene and 23.5 percent 


of stvrene were 


polvbutadiene indicated an 


investigated as to the rates of their thermal degradation. 
activation energy of 60 kilocalories per mole. 


The rates for 
The copolymer 


had very high initial rates of degradation, followed by a rapid drop, so that it was not possible 


to obtain a reliable activation energy 


Polvaecrylonitrile and a copolymer consisting of 31 percent of acrylonitrile and 69 percent 
of butadiene were investigated with regard to the nature and distribution of volatile products, 


as Wwe as to the 


rates of their thermal degradation. 


The rates were very high initially, but 


dropped rapidly so that it was not possible to determine accurately the activation energy. 
pI i | : J 


The aery 
polvbut idiene 
‘ Comparative 


rate curves 


~lvmer sample is evaporated in 35 minutes of heating, is as follows: 
| | | 


‘for NBR, 


374°C for SBR, 364° C for polystvrene, 360 


1. Introduction 


The copolymers, styrene-butadiene (GR-S) and 
aervionitrile-butadiene (nitrile rubber), currently 
designated in the technical literature as SBR and 
\BR. respectively, represent two of the more im- 
portant synthetic rubbers. This paper describes the 


results of an investigation of the relative thermal 
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stability and of the products of degradation obtained 
when these rubbers are p\ roly zed in a vacuum. A 
similar description is given of the thermal behavior 
of polyaervlonitrile and polybutadiene. 

Details of the method emploved in this investiga- 
tion have described previously i] 51.7 The 
method consists of two parts. In the first part, a 
sample weighing 25 to 50 mg is heated in a vacuum, 
first for 5 min during a heating-up period, then for 
30 min at a constant temperature ; the residue is 
weighed, and the volatile products fractionated. 
The fractions are then weighed and analyzed in the 
mass spectrometer or tested for average molecular 
weight by a microcryoscopic method. In the second 
part, a 4 5 mg sample in a platinum crucible 
suspended from a very sensitive tungsten spring 
balance [3] is heated in a vacuum, and the rate of 
loss of weight of the sample is observed at intervals. 
The activation energy of degradation is calculated 
from the initial rates by using the Arrhenius equation. 


been 


to 


2. Thermal Degradation of Polybutadiene 
and SBR 


A study of the relative thermal stability of poly- 
butadiene and the copolymer SBR, and of vields and 
chemical nature of the products resulting from 
pyrolysis of these two materials, has been reported 


previously [2,4]. Polybutadiene was found to be 
Presented before tt 2d Meetir f the American Chemical Society, New 

York, N. Y., September 7-12 

* Figures in bracket licate the literature references’at the end of this paper 


onitrile-butadiene copolymer showed a decomposition pattern similar to that of 
Here, too, the activation energy could not be deduced accurately from the 
thermal stability in terms of temperature, at which half of the 


107°C for polybutadiene, 
and 316° C for polvaerylonitrile. 


more stable than the copolymer. The ratios of small 
molecular fragments to large ones were 14:84 for poly- 
butadiene and 12:88 for SBR. 

In this investigation a study was made of rates of 
thermal degradation of these two substances by the 
loss-of-weight method, using a very sensitive tungsten 
spring balance enclosed in a vacuum [3]. Samples of 
polybutadiene and of SBR (76.5% of butadiene and 
23.5% of styrene), weighing about 5 mg, were used. 
The materials were prepared by purifying commercial 
products using a method which has been described 
previously [2]. 

Results of rate studies for polybutadiene are shown 
in figure 1. About 12 to 16 percent of the samples 
was vaporized at zero time; this vaporization took 
place during the approximately 15 min heating-up 
period. Zero time in these experiments was reckoned 
from the moment the sample had reached the operat- 
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Figure 1. Thermal degradation of polybutadiene. 
Solid lines represent percentage volatilization versus time and dashed lines 


represent 


ogarithm of percentage residue versus time, 
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ing temperature. The fact that the logarithm of 
percentage residue versus time curves deviate only 
slightly from straight lines indicates that the reaction 
involved in the thermal degradation of polybutadiene 
approaches first order. 

In figure 2 rate of volatilization is plotted versus 
amount volatilized. The apparent initial rates are 
obtained by extrapolating the straight parts of the 
rate curves to the ordinate. These rates, in percent- 
age of sample per minute, are 0.21, 0.31, 0.42, and 
0.63 for temperatures 380°, 385°, 390°, and 395° C, 
respectively. On plotting logarithm of the initial 
rates versus the inverse of absolute temperature, a 
straight line is obtained, whose slope indicates an 
activation energy of 62 keal/mole. 

Percentage volatilization of SBR versus time is 
shown plotted in figure 3 (solid lines). Here again 
there is considerable loss of weight during the heating 
up period. Plots of logarithm of percentage residue 
versus amount volatilized (dashed lines) deviate 
considerably from straight lines. The reaction 
involved in the thermal degradation of SBR appears 
to be more complicated than that for polybutadiene 
and the reaction far from first order. Rate of 
volatilization versus amount volatilized shown 
in figure 4. The rates, very high initially, drop 
rapidly with extent of volatilization. It is not 
possible to determine with any accuracy the activa- 
tion energy from these rate curves. 
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3. Thermal Degradation of Polyacrylonitrile 
and NBR 


The polvacrylonitrile * of high purity and in thy 
form of a fine white powder, had a number-averag Ay 
molecular weight of about 40,000. The NBR 
rubber was Goodrich Hyecar 1042 consisting of abou 
30 percent of acrylonitrile and 70 percent of bute {eas 


diene. The rubber was extracted for 24 hr with : 
ethanol-toluene azeotrope in order to remove im. ; 
purities. It was then kept in a vacuum at 60 to 7 pyrid 
C for 5 hr. The purified material analyzed 8.2 per} (tYos 
cent of nitrogen, which corresponds to 31 percent @ polya 
acrvionitrile. Specimens of polvacrylonitrile and! Pri 
NBR were heated prior to pyrolysis, in a vacuum at| Was ' 
110° C to a constant weight. ion 
Results of pyrolysis are shown in table 1. Frag} om 
tion V_i' represented less than 0.1 percent of the| Fre 
total volatilized and is not shown in this table. As| when 
in the case of many other polymers the ratio Vy; :Vy, color 
remained constant throughout the temperature) @"! § 
range. On the average, this ratio is 12 :SS8 in the cage) om 


of polyacrylonitrile, and 14:86 in the case of NBR. +: \ 
In the latter case, in calculating the average ratio, °, 1! 


results of experiment 1 were excluded, because vola- alter 
tilization was only 9 percent bet we 

As seen from table 1, volatilization of polyaery- Perl 
lonitrile stabilizes at about 70 to 73 percent, This and § 


is in agreement with the work of Kern and Fernow 
6|!, who found that dry distillation of this material 
Was accompanied by carbonization of the residue. A 
possible explanation of this stabilization, as given by 
McCartney [7] and Houtz [8], is the formation, along 
the chain, of ring structures involving nitrogen from 
the CN groups. NBR, which contains only 31 per 
cent of acrylonitrile, does not show any stabilizatior 
through the entire temperature range. 

Fractions V,,, from polvyacrvlonitrile and NBR 
had a tan waxlike appearance. They were hol 
readily soluble in ordinary organic solvents, including 


Phe authors are indebted to C, A. Sperat the P 
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; TapLe 1. Pyrolysis of polyacr | The main constituents in V.; are HCN, acryloni- 
| trile, and vinylacetonitrile. The composition of Vo; 
Volatile fractions» | seems to depend on the time interval elapsing 
? past on total vola- . ° ° 
1 : eeaeiiaae | 1% vite | aa between pyrolysis and analysis. This suggests the 
Experiment nur on °Re * . . . . 
tin possibility that some of the components, particularly 
v, Vos HCN and acrylonitrile, react and become non- 
es volatile on standing and, as a result, the remaining 
ACRYLONITR . . . 
volatile constituents become accentuated in the 
; analysis. This could explain why Houtz [8] re- 
, 2 2 Sh 3 : . Th, 
> 2.7 89. 1 10.9 ported that pyrolysis at 400° C yielded only a trace 
3 2 me 2.0 + of HCN. On the other hand, Hideo Nagao and co- 
165.6 86.3 13.7 workers [9] reported the evolution of a considerable 
2 43.4 sA 4 l1.¢ . oaAT > h 
fi > 52 3 a8. ¢ 4 amount of HCN when polyacrylonitrile was heated 
8 - i pe “¢ at temperatures from 200° to 350° C in an atmos- 
; 7.7 88.7 11.3 phere of air or nitrogen. Burlant and Parsons [10] 
; 7.9 87.8 12.2 . . ae 
it 7.7 0. 1 19 found ammonia as one of the decomposition products 
13 + -f- -? from pyrolysis of polyacrylonitrile. We could not 
{ 2 ‘ i = . “. : . : 
find ammonia in the volatiles. 
, SS. 1 11.9 . . . + 
, Results of mass spectrometer analysis of V2; from 
“I NBR are shown in table 3. Two of the analyses 
" were made shortly after the fractions were collected 
m4 St 15 and kept at the temperature of liquid No», a third 
; : ns ¢- 12.2 one after an interval of 3 days at room temperature, 
15.4 sis 15.2 and a fourth one after an interval of 8 days at the 
— , M4. same temperature. All the analyses were similar and 
the average results are shown in table 3. The results 
I ‘ \ -. latile at the tempera 
. } ( \ it 2 ( Dut not at 140 In i z : ’ : . 
f ott , isu found ne iry to separate this frac TABLE 3. .Mass-spectrometer analysis of Vo; from pyrolysis of 
wo pa » fac t! ter i Here it was found nitrile rubber 
\ lw ¢ 
py ridine. Average molecular weight by the micro- Formula of component A verage of 
eryoscopic method in diphenylamine was 330 +9 for { analyses 
polyacrylonitrile and 401 +3 for NBR. 
|} Fraction V_i from pyrolysis of polvacrylonitrile CH os 
was analyzed in the mass spectrometer. This frac- CoH 1.7 
e « . ry > » . CoH 3.0 
tion consisted of hydrogen. The Vey fraction C.H 29 
from NBR was too small for analysis. v= -* 
Fraction V.; from pyrolysis of polyacrylonitrile Cis 8.3 
: op 7, ne CH 3.3 
when first collected was a milky-white liquid. The | GH 2.5 
color of the liquid changed after standing, first to om te 
jtan and then to brown. Results of mass spec- ne 
° » . . . { h 
trometer analysis of this fraction are shown in table CcH 13 
2. V.; fractions that collected in experiments 1, 4, ae 4 
), 11, 13, and 14 were analyzed within 1 to 2 hours CoH 2.7 
after they were collected. The time interval C; to Cio inclusive 17.6 
between collection and analysis was 1 day for ex- rote — 
periments 2 and 7, and 1 week for experiments 3 
and Ss. 
TABLE 2 Vass-S pe trometer analysis of \ from p yrolysis of polyacrylonitrile 
rime interval bctween end of experiment and analysis 
l 2 hour 1 day 1 week 
Ex} Expt.4 Expt Expt. 11 Expt Expt.14 Expt.2 | Expt.7  Expt.3 | Expt. 8 
Mole Mole Mole Vole Mol Mole Mole Mole Mole %, | Mole ‘ 
i 4 2? 1 ois tw th IS 
\ 7] 15 44 1 7 if 
Vir { sf) 4 13 l 12 aS ( sf 59 
Py 1° 14 
Ac l 22 
Bu 9 
Proy { s f ) 19 
I 10 1M 108 1 1 100 100 100 100 100 
Exm t ‘ en in table 1 
{70183-—58 2 79 


ylonitrile and NBR 








resemble those obtained previously for butadiene and 








SBR (formerly GR-S) {2}. It is rather surprising verst 
that none of the products characteristic of polyacryl- | 5 \*} whet 
onitrile pyrolysis showed up in these analyses, al- » ple P 
though two of the analyses were made shortly after | < 
the fractions were prepared and were kept in liquid | ; - 7 
N2 until the time of analysis. It is likely that in the = 14} } oy 
presence of unsaturated fragments resulting from = 
pyrolysis of the butadiene in the copolymer, the 
hvdrogen cyanide, acrylonitrile, and vinylacetoni- | 2 | ‘ 
trile, polymerize more rapidly than in their absence, | 
and therefore do not show up in the analysis of the 
acrvlonitrile-butadiene copolymer z 4 
A comparison of thermal stabilities of the various 
polymers and copolymers discussed in this paper, 
and also of polystyrene, is shown in figure 5. The 
curves for polvacrylonitrile and NBR are based on 
data in table 1, and those for polystyrene, SBR, and : 
polybutadiene were taken from figure 1 of reference 
4}. The temperature, T,°C, at which the sample i *,* = 
loses half of its weight by vaporization in 35 min, ereeee wees guano 
can serve as a basis of comparison of thermal stability M 
of polymers. Figure 7. Rutes of volatilization of polyacrylonit 
Results of rate studies in the tungsten spring bal- 
ance are shown for polyacrylonitrile in figures 6 and | peratures, a considerable initiation period is 9}! 
7. In figure 6, percentage of volatilization is plotted | served. Beyond this period volatilization proceed 
versus time from start of experiment. At lower tem- | rapidly, but soon drops to a very small value, ing. 
cating stabilization of the residue. At higher tem. 
peratures the initiation period terminates during th 
. 15-min heating-up period and the curves show stabj. 
ws * ization of the polymer soon after the sample reaches 
ow | the operating temperature. These facts can be ob} 
‘did ‘ served more easily in figure 7, where rate of volatij 
' zation in percent of the original sample per minut! 
| is plotted against percentage volatilization. It is ne 
possible to estimate with any degree of accuracy th! 
; initial rates or the activation energy based on a 
of the rates shown in figure 7. 
a —7 ae oe —. Results of rate studies for NBR are shown plotted 
EMPERATURE , DEGREES CENTIGRADE in figure 8, where percentage volatilization is plotte 
FPigure 5 Relative thermal stability of polymers and 
copolymers 
Cu poly “BR “ bu ‘ i 1 j 
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| versus time from start of experiment, and in figure 9, 
«here rate of volatilization in percent of original sam- 
» ple per minute is plotted versus percentage volatili- 
f mation. The two sets of curves resemble closely those 
} for SBR shown above (figs. 3 and 4), and no definite 
conclusion as to initial rates and activation energy 
} ean be based on these curves. 
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4. Discussion 


The main constituents of the light volatile fraction 
V,; obtained in the pyrolysis of polvacrylonitrile, are 
tvdrogen evanide, acrylonitrile, and vinylacetoni- 
trle. The hydrogen cyanide most likely forms 
by splitting-off the chain in the same manner 
as CH,COOH in the pvrolvsis of polyvinyl acetate 
11], HCL in the pyrolysis of polyvinyl chloride [12], 


ind HF in the pyrolysis of hydrofluoroethylene 
polymers [13]. 
H H H i H 
( ( ( ( If ( ( 
Hy Aa H C/H H 
( ( ( > ( ( HHICN 
CN CN CN CN CN l 


Formation of monomer can be visualized as follows. 
} \ random SCISSION of n (' & bond in the chain re- 
} sults In two free radical ends. 
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H H He} I II H 
( ( ( : ( ( ( 
H/H\ H/H\VA/H\IAY/EHY/HR\H/aH 
~ ( ( '¢ ( Cc > 
CN CN CN CN CN CN 
H He} H Hq } H H 
( ( ; Cc: Cc i Cc Cc 
H/UNH/uU W/u H/UHN IH Y/R BY AS, 
( ( ( “Cc °C Cc 
CN CN CN CN rT N CN 
iH H H i 
( C- C Cc 
HyUNHYH iH H H H H/UN WSU 
( ( - Cc C C+-C C 
i iH | 
CN CN CN CN CN CN (2) 


This is followed by an unzipping 
monomers, as shown in eq (2). 

Formation of vinylacetonitrile could be explained 
if we assume a head-to-head arrangement at some 
positions in the chain, as shown in eq (3). 


reaction vielding 


1 CN 

1 i ( 2 H H 

( ( WN 17°C c C 
HI TL“ ®/H\W/H rT 
( ( { q> (  * Gc > 

HI CN CN CN 
CN CN CN 
2 free radicals + CHy=CH—CH;—-CN (3) 

This formation involves 2 (—C' secissions in the 


chain, followed by a hydrogen transfer from position 
1 to position 2. 

The pyrolysis curves for polyacrylonitrile (fig. 5) 
as well as the rate curves (figs. 6 and 7) show stabili- 
zation during pyrolysis. This is characteristic of 
polymers in which the loss of side groups from 
adjacent carbons of the chain results in conjugated 
double bonds in the chain, since these double bonds 
tend to enhance the thermal stability of the chain. 
Polyvinyl chloride [12] and polyvinylidene fluoride 
[13] may be cited as examples of a behavior similar 
to that of polyacrylonitrile. 
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' On the Motion of Two Cylinders in an Ideal Fluid* 


Lloyd H. Carpenter 


The complex potential of two cylinders moving in an infinite liquid is determined by the 
method of image doublets, and the solution is expressed as an infinite series in rectangular 
coordinates. Approximate solutions in finite form are given for various cases. A method for 
generalizing the solution for the case of more than two cylinders is indicated. Applications 
to the flow induced by a cylinder moving in the presence of plane boundaries are given and 
the stream lines are illustrated in certain cases. 


1. Introduction | y 


The problem of the motion of two cylinders in any | 
manner with their axes always parallel was appar- | 
ently first solved by Hicks." The velocity potentials | 
are found in general as definite integrals, which, | 
when the cylinders move as a rigid body, are ex- | 
pressed in finite form as elliptic functions of bipolar 


Z=xtly 


"SS 








coordinates. Greenhill? has expressed the solution | ; >x 
as an infinite series in bipolar coordinates. 
In a problem in wave motion the need arose for a | 
solution of the potential flow for a cylinder moving | 
parallel to a plane boundary. The problem is the | 
| same as that of two parallel cvlinders of equal radii | 
moving with equal velocities normal to the line of 
centers. Thus the flow iS a special case of the classic FiGuRE 1. Notation diagram for motion of two cylinders in an 
results obtained by Hicks and Greenhill. However, infinite liquid. 


it was necessary to have the potential function 
represented in rectangular coordinates to determine 
the additional term representing the wave motion. 
The bipolar coordinates used by Hicks and Greenhill 
are defined geometrically, and the transformation to 
rectangular coordinates does not admit a convenient 
representation. Therefore, the solution was _ ob- 
tained by the method of image doublets and is here 1of0'\?  e-ta ‘ 
generalized to the case of two evlinders moving in any Ww, { b?( f ) Pay | f.) (2) 
manner with their axes always parallel. The present ie 
method may also be applied to the motion of more 


where 2 is the complex coordinate position. The 
normal velocities introduced on C’ by w are cancelled 
by the image doublet whose axis is in the direction 
Ta 


where 


than two cylinders which was not done in the bipolar » FP 
system with either Hicks’ or Greenhill’s analysis. " f° (3) 
Only ideal fluids without vorticity are considered 
in this study. The sum w+; is the well known expression for the 
potential of a doublet in the presence of a cylinder 
2. Solution for Two Cylinders (see Milne-Thomson *). Similarly, the normal ve- 


ae allel evlinders C and (’ of radi locities introduced on C by w, are cancelled by the 
onsk er tWO para el cy inh ers ari oO re il second image doublet 
band 6’ momentarily at a distance f apart moving 
normal to their axes in a liquid at rest at infinity. b’\2/ b \2 ote 
Nn . . ; . 4 rhe ) ) é 
Choosing the coordinate system as in figure 1, the w.—l b? ( -) ( if ) (4) 
evlinder velocities may be expressed as Ue'* and 7 : 
U’e'’, respectively. 

The complex potential uy Py iWo due to C, in where 


the absence of ©’, is that of a doublet whose axis is in h2 
the direction a jo=>" (5) 
I—Ji 
Ub? l 
wr - ’ (i) ‘ . . . a . 
. Continuing, an infinite sequence of image doublets 
. ; . : - 

———_——_——. of decreasing strength is obtained. The general 

i his problem was encountert n connection with a wave investigation spon : a ° f 2 = 
wred by the Office of Nava Researel ind the present account is published in expression Or W, 1S 
rder to make more convenient formulas available 

j W. M. Hicks, On the motion of two cylinders in a fluid, Quart. J. Pure and 

Appl. Math. 16, 113 (1879 
,, A. G. Greenhill, Functional images in Cartesians, Quart. J. Pure and Appl L. M. Milne-Thomson, Theoretical Hydrodynamics, Article 8.82, p. 221 
Matt 18, 231 ISS2 MacMillan & Co., Ltd., London, 1938). 
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=_— ) 3. Approximate Solutions 
B;. - ; for k even | 
Wy 4 , > 6 The series | 12) converges rapidly for small valusf W 
eta : of bb’/f so that for this case one or two lerms of thel nfin 
B —(f—f \for k an series gives a close approximate solution ana 
ie v- 
} if 
h’ , ; , f bb \] ne 
_— (0 wW=W,4+- Ww w Ww or. small. y, OH 
where B—l'b?, B —( b? ( ; ) ’ _ %) eoor 
. . ~ \ the | 
hh? For larger values of bh’ /f a more accurate solution! loul 
3 ae 3 9 7) may be obtained by approximating the remaining! the | 
7 oS a terms. For this purpose let *| befo 
and . =" 
n=l|im f, (k odd lim f,” (k even ers 
D > . > man 
}; > l4 
for k even | n’=lim f, (k even) =lim f,” (k odd) loul 
} . 
. > > start 
/ 4 hy r ia i , (vai 
for I: odd rhe conditions f th 
LJ f th 
j n )n h 
it being understood f;=0. ; 
/ n)” h lf 
Similarly, starting with the potential of C’ in the If 
absence of ¢ F uid 7» * b, 7 h ia be used to det rine y andy etw 
Now JOF / greater than some number 7 replac solut 
7 j hy ” fol / evel by n lol / odd. and replac mov 
w l bh’)? ;’ ‘ f," by n for & even, by n’ for k odd giving » she 
hp iwa 
b ‘ rb " - | TI 
° , , ( ) . / ou j 
w h ( ;) f? 4 / n f n : ~ > move 
erel 
iLtit and omy] 
, ( ‘ , / Dass’ 1 { 
D for | evel ul Du 9 l\ 
, e / / 
” 10) ; : _ , i 
il . for & sufficiently laren Thus the solution from y 
re] j fo} j odd would be 
J 
\) 
where 4 8 
lh u u u Th / u erm 
bud 3 
, ; , j Lei 
p [ h 5 / h ( ) 
Me CIUuLSE 
hh’ : 
D p i] , 
t{—f {—f SS‘, , y, 
. cy l—, 
and 
‘a h An accurate solution for any value ol bb’ f can 
for k even obtained from eq 19) by an appropriate choice ol 
/ ) lla lt Wiis found It) the computatiol : for the examples t 
h? j follow that for bd’ / 6, | gave an accural 
/ / lor k odd J result, and for the extreme case bb’ f , an accural The s 
solution was obtained for n » lato 
and again f,’=0 Again when bb’) is sufficiently small » and at foyy, 
Thus the complex potential uw ) YW for the negligible so that eq (19 reduces to / 
motion of two evlinders is Th 
a ail | er | 
Pi u uw w u u hh’ ” 
‘ ) - “qua 
u S u u 12 | ( , ) | 
. i , : 7, 
where the w, are obtained from eq (1) through (7 where wp, Ww), Wy’, and aw,’ are given by eq (1), © 
and the w,’ from eq (8) through (11 8), and (9) after setting f/f, and /,’ equal to zero 


84 








4. More Than Two Cylinders 


nal Valet When more than two cylinders are moving in an 
IMs of the jnfinite liquid the method of solution is completel) 
analogous to that given in section 2. ; 

} If n evlinders are moving in a liquid choose any 

one Cp as a reference and introduce an appropriate 

: 3} coordinate svstem. The complex potential of C, in 

the other evlinders will be that of a 

The image of wp. Is taken in each of 

| evlinders in the same manner as 

these »—1 tmage doublets will 

eg and in each of the other eylin- 

lers as well The continued in this 

manner with each doublet introducing n—1 image 

4) Joublets at each step. The same process is applied 

! starting with each of the hi | remaining evlinders. 

vain the potential of the system will be the sum 

f the potentials of all the doublets. An application 
{this method will be given in section 6. 


the absence ol 
} ] . 
le Solution! loublet wp, 
ithe remaming 7 
hefore Each ol 
ihave an image 1th ( 


remaining! 


process Is 


5. Cylinder Near a Plane Boundary 


If b 5. U’ . a’ Tr a the plane equidistant 
hetween the evlinder axes is a streamline and the 
flow induced by a evlinder 
moving in the presence ofa plane boundary Figure 
? shows the streamlines for a evlinder moving toward 


mn andy 
vv replac 
dl replace 


solution represents the 


iwall, a=0 
The equations can be simplified when the evlinder 
}moves parallel to the wall, a=2/2, with its cireum- 
erence touching the wall, bf 1/2. In this case the 
omplete solution Is 


irom . 7 
VP k+l hb fh lJ 
3 )An approximate solution is obtained by taking » 
i terms of the series and in the remaining terms re- 
p lacing keh | and (4 2 I | bry | Thus 
; ‘4 | ‘ 
lbs | r | 
fd (k+-1 4 h 2: ki? \ 
y | / ! as : | 
2 L k+1 6 & Lj 
» fe 
f° can 7, a or 2 
oice ol 14 0 fad (h 1) 23 
amples /, | 
Linh pie L J 


accural 


accu The streamlines are shown in figure Ba. The compu- 
‘The streamlines for 


lations Were made with n—3. 


nd 9 a foure 3b were obtained by superposing a velocity 
(on the evlinder and the liquid. 
The streamlines for liquid flowing around a eylin- 
ier parallel to a plane wall are shown in figure 4 
2) Equation (19) was used with n—1, a=2/2, b/f=1/4. 
i 
(1), @ 
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Ky URE 2 Streamlines for culinder moving loward a plane 
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FrGcurE 3b Streamlines for current flowing with general 
ocity U past a fired cylinder against a plane boundary. 
Obtained from figure 3a by impressing a velocity U on the fluid and the 

{ 1, ? l 
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FIGURE 4 Streamlines for current flowing with general velocity 


I past a fixed « wlinder neara plane ho in lar / 








{ 1, b=1 
6. Cylinder Near Two Plane Boundaries | @qqeqoyee ee 
> , ' - ; | Fieure 6 Streamlines for a cylinder mor iq toward one oft i — 

If four cylinders of equal radii have svmmetric perpendicular plane hae 2 eas aa’ 
positions and velocities with respect to two perpen- U=1,! 
dicular planes as indicated in figure 5, then each = 
plane is a streamline and the solution represents the : , — — 
flow induced by a evlinder moving in the presence 7. Cylinders Moving as a Rigid Body —_ 
of two perpendicular plane boundaries. Consider- if ths welactty. Tet. j : 
ing the case f=4a, g=4a, a=0, an approximate solu- t the velocity, Ue’, is the same for each cylinde,) — 
tion is obtained by the method of section 4 stopping | [N@Y Move as a rigid body. ‘The previous solution} —_ 


for a cylinder moving parallel to a plane boundary} 
is one example of this case. Another example 5} 


given in figure 7 where b’=b=f/2, a’=a=0. Th! — 
solution simplifies to | 


after the first set of images. The streamlines are 
shown in figure 6 





. ! Pp , ‘ 1)" I } 
w=Ub 3s — - gal Fiat 
7 (k+1)?< 2 ke > k+2 24 
6 k4+1 Ob kt > 
’ C2 The approximate solution obtained by the metho 
; of section 5 is 
2 
— _—_ a = 2 a s - < _ { 2 r ‘5 
F wutb >> — ae 
: <I (k+1)? 2 2 k > k+2 
ye \ 6 k+1 6b k+l] 
{ma 
| ym, \ ) k ‘ 
| \ ™ / 12 0 (k +-])* 25 
“3 \VA 4 } | 
j ) . 
¥ J 
U > 


The computations were made with n=3. 
( 5 lat 1 diaaran / mol ; iacy rele ner ro : ~ : 4 } 
Ficure 5. Notation diagram for motion o ; ato Another example is given in figures Sa and § 
per pe ndicular plane houn ! ‘ h’ } . , vrN ’ ° 1 
where 6’=b=f/4, a’ =a=0. The solution was 0b- 

tained using eq (19) with n=1. 


Che dotted lines indicate the streamlines that can be cor lered as solid bound 
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Streamlines for current flowing with general velocity 
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Streamlines for current flowing with general 
velocity U past fixed cylinders. 
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Research Paper 2890 


Properties of Sodium Titanium Silicate Glasses’ 
Edgar H. Hamilton and Given W. Cleek 


The glass-forming area of the Na,O-TiQo-SiO2 system was surveyed. 
th refractive indices from 1.5184 to 1.8005, nu values from 51.5 to 23.2, and densi- 
The glasses have higher dispersions and have considerably lower 
the PbO glasses with comparable refractive indices. 


2 42 to 3.00. 


densities than 
1. Introduction 


In order to give the optical designer greater free- 
dom in his choice of glasses for use in optical devices 
there is a need for glasses that cover a considerable 
range of refractive indices, and for each index value, 
arange of dispersions. Most commercially available 
optical glasses having high refractive indices and 
high dispersions contain PbO in varying amounts. 
In general the refractive index and the dispersion 
increase With increasing PbO content. However, in 
PhO glasses there is a limiting value of dispersion 
fora given value of refractive index. 

It has long been known that small additions of 
TiO, to a glass causes a relatively large increase in 
refractive index and dispersion. One of the early 
references to the use ol ho in lass is that of Har- 
IS70. He is reported to have produced 
“two ‘nearly flawless’ 3-inceh disks of titanium glass 
and two of borate glass” [1].- The well-known work 
of Schott and Winkleman [2], in association with 
Abbe, the optical designer, was the first S\ stematic 
attempt to determine the effect of chemical composi- 
tion on the physical properties of glass. The imme- 
diate object of their work was the development of 
glasses having desirable optical properties. ‘TiO, was 
one of the 28 oxides for which they studied the effects 
ol composition on the phvsical properties of glass. 

Colbert [3] has compared the effect of TiO, and 
PhO on the refractive indices of glasses. He added 
TiO. to a base glass of the composition Na0,351¢ ), 
and PbO to a base glass of the same composition 
and found that in each series the refractive indices 
nereased sharply 

Three in the ternary system, 
\a,0-TiO,-SiO, were reported by Turnbull and Law- 
1]. However, they did not survey the glass- 
lorming area in the system. More recently Varguine 
5} has reported on the glass-forming area of the 
\a,0-TiO,-SiO, system, but he gives data on only 
one series of glasses in the ternary system. 

Thus, although studies have been made of glasses 
containing TiO,, it was deemed desirable to make a 
systematic survey of the entire glass-forming region 
of the Na,O-TiO,-SiO, svstem, because it appeared 
1 this ternary system could be 
ised for the development of optical glasses having 


court about 


of glasses 


Series 


ence 


that compositions | 


high refractive indices and high dispersions. Some 
Dinilietenatensann 

Che work d bn ! \ ed by the Bureau of Ordnance. 
’ ment « thy N\ 

Figures in | ket fe the end of this paper 


Glasses were 


of the properties of the glasses made in this survey 
are reported. The phase equilibrium diagram of the 
ternary system has not been determined, but infor- 
mation on the binary sides of the ternary system is 
available from previous work [6, 7, 8}. 


2. Experimental Procedure 


The glasses were made in 500-¢ melts from batch 
materials of sufficient purity to satisfy the require- 
ments for optical glass. The melts were made in 
platinum crucibles and stirred with a platinum—10 
percent-rhodium stirrer. The melts were cast in the 
form of blocks approximately 3 by 3 by % in. The 
details of the melting procedure are given in previous 
publications [9, 10]. The sag point [11] of each glass 
was used to establish an annealing temperature from 
which the glasses were cooled at 245° C per hr to 
350° C 

The liquidus temperature of each glass was de- 
termined by a temperature gradient method [12]. 

Refractive indices were determined on polished 
samples of the glasses in the form of 60° prisms for 
the C, D, F, and G’ lines by the NBS Optical In- 
strument Section. 

The densities of the glasses were determined by 
the method deseribed by Glaze, Young, and Finn 
[1:3]. 

The glasses were not analyzed. The compositions 
in table 1 are all caleulated from batch. 


3. Results 


3.1. Glass-Forming Area of the Na,O-TiO,-SiO, 
System 


The compositions of the melts are given in table 1 
and are plotted in the triangular diagram of figure 1, 
Binary sodium silicate melts were made containing 
from 15 to 40 mole percent of Na,O. The glass 
containing 15 mole percent of Na,O devitrified during 
annealing and the other binary glasses developed a 
hazy appearance. The haziness of the annealed 
glasses increased with increasing SiO, content. Due 
to their hygroscopic nature none of the binary 
classes could be polished in a water medium. This 
was also true of glasses containing 5 mole percent 
of TiO, and 25 or more mole percent of Na,O. All 
the other ternary glasses containing 35 or more mole 
percent of Na,O could be polished in a water medium 
but did not majntain a polished surface for more than 
a few weeks time 
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E1973 25 
E1974 30 
E1977 33.3 
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E1976 19 
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E1905 a0) 5 
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EF 1968 35 5 
E1069 10) F 
E1970 i 5 
E1922 10 10 
E1999 12 1) 
E1960 15 10 
E2002 17.5 10 
E1961 20 10 
£2003 22. 5 10 
E1907 25 10 
E2000 97.5 10 
E1962 30 10 
E2001 32.5 10 
E1908 35 10 
E1963 th 10) 
F.1807 45 10 
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E1925 10 l 
E1987 12 1 
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See footnote at end of table. 
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TABLE 1. Ternary Na,O-TiO,-SiO, compositions— Continued 


Melt NwoO rio SiO nD , p Lquidus Primary Remarks 
temperature) phase ! 
Vole Vole Mole ( 
E1953 20 i 1 k Devitrified in mold. 
h } E2027 ») i ; 2. O84 
E1954 2 ' 35 1. 8005 23. 2 2. 995 1060 I Devitrified streaks in glass. 
F195 4) j 30 1.7874 23.4 2 GSS 112 b Stones, 
E1956 5 i 20 Ys0) F Devitrified in mold. 
E1958 , , du) 1102 I Devitrified in mold. 
} E1957 il . 2 1040 k Do. 
{ 
E1883 2 4 25 Devitrified in mold. 
l. 
Primary pha \=SiOe; B= Na SipQs; C= NapSiOg; D= Unknown; E= NajoO-TiO:-SiOg: F Unknown; F)= Unknown. 
MPOUN 4oor 
R ” 5 MOLE No 
ry ; 7 : 
. Ln RG r x F 
AY F F 30 MOLE % Na, 
} W — 
. Ww a \ y “= 
‘ WwW } . A. 
| MOLE % No x 
n * s r Se 
uw ‘eo wo 
x 8 —_ 
Je ws 
. . } 5 MOL & . 25 MOLE % No.0 
ration eee . ‘ x ; os No, MOLE © NOQV 
. ° ; " 
« 2 Fs aw 
° a - 
° x Fad F i we 
———/ 
= ro £ od 
, ——s 
> 0 20 30 40 5 y 20 30 40 50 
MOLE PERCENT TiO, 
Figurt 2. Liquidus temperatures of glasses in five series of 
Na.Q-TiO,-SiQ» glasses. 
1 Primary phases present at the liquidus are: A, SiOe B, NagO-28iO," C, Na,O- 
SiOy E, NayO-TiO:-SiOe D, F; and F) were not identified. 
} 
the crystals appearing at the liquidus and from 
) Ficere 1. Glass-forming area of the Na,O-TiO.-SiO, system | 20 to 50 degrees below the liquidus, are too small 
| Primary fields are: A, SiO» B, NasO-28i0e C, NaO-Si0e E, NaO-Tids-sior | to permit the determination of their properties by 
D, Fy and Fo, unidentific 


microscopic means. At lower temperatures phase 
All of the ternary glasses containing 5 mole | D is observed, but it is believed that the submicro- 
percent and some containing 10 and 15 percent of | scopic crystals are not phase D. Adjacent to the 

| TiO, were practically colorless except when viewed | phase F; field crystals with the same appearance and 
edgewise. The remainder of the ternary glasses | properties with the exception of refractive indices 


were either amber or yellow in color. As the | as phase F,, were observed. 

Na,O/SiO, ratio is increased, the amount of color The properties of the crystals from phases D, E, 
in the classes decreases rapidly, so that glasses are and F, are: 

obtained whose color compares favorably with PbO Phase D. Composition unknown. Thick rods 


glasses having comparable optical properties. The with high birefringence, am x about 1.70 and Main 
composition ranges of the hygroscopic, amber, and | about 1.60. 


vellow glasses are indicated in figure 1. Phase ©. Composition believed to be Na,O-TiO,- 
SiO,. Rod-shaped crystals with low birefringence, 
3.2. Liquidus Temperatures parallel extinction, negative elongation. Mmax about 

The liquidus temperatures of m: ny - the glasses | 1.77 and Nain about 1.73. 
are given in table 1 and plotted in figure 2. The Phase F,.. Composition unknown. Needle-shaped 
letters under the curves indicate the he imary phases | crystals with moderate birefringence, parallel ex- 
at the liquidus. The letters, A, B, C, and E stand for | tinction, positive elongation, Mx about 2.00, and 


SiO,, Na,O-.2SiO,, Na.O-SiO., and Na.O-TiO.-SiQ,, | Main below 2.00. 

respectively. The compositions of the primary 

phases designated D, F,, and F, were not determined. 3.3. Refractive Indices and Dispersions 
The approximate composition ranges of the primary 


fields are indicated, insofar as the ‘vy are known, in Glasses were obtained with refractive indices from 





figure 1. A region designated F, is located adjacent | 1.9184 to 1.8005 and nu values * from 51.5 to 23.2, 
to the D and F primary fields. The approximate | — a 

) boundaries of this field, if it is a primary field, are italien i a 

|; hot known. In the F Geld adjacent to the D field | ; — ~— 
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The refractive indices of the glasses are plotted in 
figure 3. The lines connect glasses containing equal 
moles of Tif ». No simple relation is readily evident 
between refractive indices and in Na.O 
content of the glasses. 

The nu values of the glasses are plotted in figure 4. 
As in the the refractive index curves of 
figure 3, the nu value curves change slope and direc- 
tion as the TiO, content of the glasses is increased. 


Increase 


case of 


1.82 
1.80 ~ 
i. 
“40M 
1.78 
a 
1.76 “ 
~_ 
~. 
1.74 } i M 
| 
| e?*-s 
1.72 } , 
| —_ 
oO | > M 
c |.70} 
= | 
x o*——6. 
wy 1.68 } Pints ‘ 
= 1.66} 
| 
uJ | -_ 
| 4 » se 
> 6 x i M 
i ¥ 
© 1.62} 
<a 
os | as 
u_ !.60} a> -_ s+ « V A 
WW | - soll 
x | 
1.58 } 
- ” 
56 | pe 
' . 
| aa 
54 | 
> mM 
52 | ° 
ee 
5 . xX c 
M RCE 
V " 
Figure 3. Refractive ind of the Na,QO-TiO.-SiO, glasses 
(riAasse which contained striae are identified bw * nd hvgroscopic%¢ ' 
by 
= 
& 2 un, 
~ | 
| 
4a. 
2 > > ——.@— 
4C 
j e -—- - = 
Ss | 
— | 
aq 350} 
> — ° 
rn a aA 
2 
e ° 
30.C —we 
> 2 ‘ 
~ as . 
- 
- ~f" 
* - > 
— » 
eo —_* 
= ah 
aA 
. 5 > fa 
M E ERC 
Fig ! Vu ! ( j Ni®O-TIOSIO F 
(ria \8 t j 


92 





In figure 5, the optical properties of ty pical opticg) 
flint glasses [14] and some of the Ni.O-TiOeg 
glasses are compared. All of the curves for the eens, 
imental glasses seem to parallel each Other, The 
glasses containing 25 and 30 mole percent of Nag 
all plot on the same line. Ty pical optical flint gla . 
with refractive indices above 1.605 have lower dis} 
persions (higher nu values) than the Na,O-Ti(,. 
SiO, glasses with corresponding refractive indices — 
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3.4. Densities 
The densities of the classes containing PiO, are 


considerably lower than for PbO classes with the 
same indices of refraction The density of the Tid, 
glasses range from 2.42 to 3.00 while the density of 
the PbO classes cover the range of 3.0 to 5.0. For 
this reason the less dense TiO classes may be of 
interest in applications where weight is a factor of 
prime importance. 

The densities of the classes are plotted in figure 6. 
There is a definite increase in density from series @] 
series as TiQ, is increased. As Na.O is increased i 
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Wavelengths From Thorium-Halide Lamps 
William F. Meggers and Robert W. Stanley* 


The present system of international secondary standards of wavelength for spectro- 
scopic measurements is based on interferometric determinations of wavelengths emitted at 
atmospheric pressure by an electric are between iron electrodes. Because of the poor quality 
and uneven distribution of these iron standards they are not suitable for accurate measure- 
ment of wavelengths in the spectra of heavier elements, most of which are more complex 
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and consist of much sharper lines than the standards. 
when excited by 


quantity of a thorium halide, 


Quartz-tube lamps containing a small 


microwaves, emit thousands of uniformly 


sharp and evenly distributed lines whose wavelengths, or positions in a spectrum, can be de- 


termined with about one-tenth the 


oo 


interferometers with plate separations of 25, 


grating spectrograph in making these measurements. 
and from 
value of 27 classified thorium lines is tested by means of the combi- 
average error is less than 1 part in 20 million. 


3288.7356 to 6991.5839 
The accuracy in relative 


nation principle, which indicates that the 


A in vacuum 


l. Introduction 


During the past half-century the wavelengths of 
hundreds of radiations from the iron are at atmos- 
pressure have been measured with inter- 
ferometers to serve international standards of 
wavelength for spectroscopic measurements. The 
history and results of this activity were summarized 
in Wavrele ngths From Tron-Halide Lamps [1] I in 
which it was shown that the accuracy of measure- 
ment could be increased by a factor of 2 or 3 if the 
iron lines were emitted by simple lamps containing 
a trace of iron halide excited, at low pressure and 
moderate temperature, by microwaves. But the 
iron are in air is still the main source of interna- 
tionally adopted standard wavelengths. 

The actual application of the adopted iron stand- 
ards to wavelength measurements in various spectra 
s handicapped by two serious defects in the iron-are 


pheric 


as 


spectrum. First, because of the high temperature 
6,300° KX.) of the iron are and the relatively small 
tomic number (Z=—26) and mass (A=56) of the 


ron atom its spectral lines are excessively broadened 
m account of Doppler-Fizeau effect. Moreover, an 
deetrie are between iron electrodes at atmospheric 
pressure produces lines of various kinds; strongly 
self-reve rsed lines involvi Ing low ene rey leve ‘ls, diffuse 
lines involving high energy levels, and unstable lines 
f variable wavelength due to “pole effect.’ Second, 
the distribution of strong or similar lines in the are 
spectrum of iron is such that no uniform spacing of 
wavelength standards is possible, and unduly large 
zaps between usable standards are unavoidable in 
ertain spectral regions. Consequently many spec- 
inoscopists have expressed dissatisfaction with the 
iternational standards of wavelength emitted by 
the iron are. Thus, Harrison [2] states that the 
vavelength precision of measurements in the MIT 
Wavelength Tables has been limited primarily by 


LT 
* Present addre Purdue University, Lafavette, Ind 
Figures in brackets indieate the literature references at the end of this paper 
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error of locating iron-are lines. 
vacuum wavelengths emitted by a thorium-iodide lamp have been measured relative to 
5462.2705 and 4047.7144 A emitted by a similar lamp containing mercury-198. 
10, or 50 millimeters were used with a stigmatic 


Preliminary values of 
Fabry-Perot 


The thorium wavelengths range from 
3287.7885 to 6989.6562 A in standard air. 


“the insufficiency of adequate wavelength standards 
in some spectral regions, particularly of standard 
lines of suitable intensities.”” Likewise, Bovey [3] 
in giving a preliminary list of 6-figure wavelengths of 
plutonium excited in a furnace ‘complains that “A 
serious difficulty for accurate measurement has been 
the width of the iron-are lines, many of which are 
several times wider than the emission lines from a 
furnace.”’ Also Fred and Tomkins [4] in describing 
the first two spectra of americium to 7 figures, an- 
nounce that ‘‘Most of the uncertainty in the wave- 
lengths is due to the poor quality of the standards, 
since the agreement in the hfs intervals shows that 
the line positions can be measured to another decimal 
place.” 

Logically the iron-are standards should have been 
used only for measuring atmospheric are and spark 
spectra of metals lighter than iron, because, other 
things equal, the heavier atoms emit sharper lines 
which cannot be measured accurately relative to 
standards that are broader and badly spaced. When 
the are spectrum of cerium (Z=58, A= 140) is photo- 
graphed beside that of the standard iron are it is 
readily seen that all cerium lines are narrower than 
iron lines, and a comparison with thorium (Z=90, 
A 2) or with uranium (Z=92, A=238) is even 
more embarrassing for iron. This was observed by 
Harrison |5|, who remarked that “In the extensive 
series of wavelength measurements made with the 
automatic comparator, which resulted in the publica- 
tion of the Massachusetts Institute of Technology Wave- 
length Tables, the most serious limitation on precision 
was found to arise from the lack of satisfactory stand- 
ard lines. Measurements on lines in such complex 
spectra as cerium and thorium, for example, showed 
much better internal consistency than did those on 
the iron standards themselves.”’ To assist in remedy- 
ing this situation, Harrison [5] cleverly contrived a 
complicated machine for the rapid direct determina- 
tion of wavelengths from Fabry-Perot interferometet 


patterns. He called this machine the WINMAC 








(Wavelength interferometric measurement and com- | 


putation), and in 1950 he reported [6] that ‘‘the 
WINMAC is now being readied for what is hoped 
will be an extensive series of wavelength determina- 
tions of intermediate wavelength standards in cerium, 
thorium, and other elements.”” Unfortunately this 
hope did not materialize. However, this is not seri- 
ously regretted now because in the meantime great 
improvements over the traditional ares and sparks 
at atmospheric pressure have been made in spectro- 
scopic light sources. 

The most simple, convenient, efficient, and eco- 
nomical type of spectroscopic light source is an evac- 
uated tube of fused quartz enclosing a trace of gas, 
metal vapor, or volatile metallic compound excited 
by microwaves generated by vacuum-tube oscillators 
or by magnetrons. The prototype was the low pres- 
sure, water-cooled Hg'’ lamp deseribed by Meggers 
and Westfall [7], who also discovered that “clean-up” 
in such discharges varied inversely with frequency. 
In most laboratories magnetrons producing 2,450 Me, 
more, soon replaced vacuum-tube oscillators pro- 
ducing 100 300 Me for the excitation of these 
lamps. The preparation of similar lamps with other 
volatile metals was described by Zelikoff, Wvekoff, 
Aschenbrand, and Loomis [8], while Corliss, Bozman, 
and Westfall [9] demonstrated that such lamps can be 
made to emit the spectra of all metals that form vola- 
tile compounds. <A precise procedure for the prepa- 
ration of electrodeless discharge tubes containing very 
pure rare earths and highly radioactive elements has 
been described by Tomkins and Fred {10}. Conrpared 
with conventional ares and sparks, these lamps re- 
quire a minimum amount of sample, which is enclosed 
and conserved. Also, they operate at relatively low 
temperature and pressure, and consequently emit 
spectral lines of greatly reduced width usually free of 
self-reversal. When operated with continuous wave 
power and moderate gas or vapor these 
lamps strongly favor the spectra of neutral atoms, 
but pulsed discharges and/or reduced pressures en- 
hance the second and third spectra. Lamps viewed 
end-on when excited by a 100-watt magnetron emit 
spectra comparable in intensity to a 10-ampere 220- 
volt d-e are, but with line widths much less than half 
Bright. sharp lines ure especially needed 


to 


pressure, 


as great. 


for the resolution of complex Zeeman patterns, or for 


the measurement of isotopic-spectroscopic effects 
The development of these electrodeless lamps has 


thus greatly increased the range and precision of 


spectroscopic observations 

Now the accurate description and quantum inter- 
pretation of the spectra of rare-earth types of atoms 
and ions is still in a very unsatisfactory state, because 
of their extreme complexity and because pure samples 
have not been available in the past. The latter 
difficulty has recently been overcome either by 
efficient separation or by artificial preparation, but 
extreme complexity of the spectra remains as an 
It is primarily for 
this reason that wavelengths must be measured with 
greater accuracy, and for this purpose reference 
standards of greater precision and better distribution 
than those from an iron are are absolutely necessary. 


obstacle to successful analysis. 





Fortunately the simple electrodeless lamps referred 
to provide the means for obtaining a greatly im. 
proved system of standard wavelengths [11]. Lamps 
containing a trace of a thorium halide excited by 
microwaves emit approximately 20,000 lines of yj. 
formly sharp character and even distribution be. 
tween 2000 and 12000 A. Since these lines are! 
characteristic of heavy atoms and ions of even mage 
(7=90, A=232) they have relatively small Doppler. 
Fizeau widths, and are entirely free from hyperfine 
structure or isotope shifts. By imaging a thoriun. 
halide lamp inside a similar one containing Hg" (o} 
vice versa) both sources can be caused simultaneously 
to illuminate a Fabry-Perot interferometer and 
stigmatic spectrograph. From the interference pat. 
terns the wavelengths of thorium lines can then be 
determined to 8 figures relative to the previously 
determined values of the mercury lines. The pyp. 
pose of this paper is to present preliminary results of 
such wavelength measurements for 222 radiations 
from thorium, covering slightly more than one octaye 
of spectrum. 





2. Light Sources 


The thorium source used in this investigation Was 
an electrodeless discharge tube containing a fey 
milligrams of thorium iodide and pure helium gasat! 
7] pressure of 5 nim of He or less The tube Was pre- 
pared by C. H. Corliss in a manner previoush 
deseribed 19] The tube blank is made of V year | 
glass; it has a bore of approximately 5 mm anda 
length of 11.3 em. The discharge was excited at a} 
frequency of 2,450 Me obtained from a commercial] 
diathermy venerator. } 

The thortum tube used throughout the investige 
tion was chosen from among several that were mad 
at the same time. After a warmup period of 30 tof 
the tube would glow very brightly for approx 
mately > min At the end of this time the thorium| 
spectrum would weaken, and the discharge would be 
confined to a narrow filament near the center of the 
tube, indicating that the pressure in the tube had 
increased. The period of high intensity mentioned 
above was accompanied by a gradual movement d 
the thorium iodide deposit from the top to the bot- 
tom of the tube. As soon as the thorium spectrum 
began to weaken, the discharge tube was removed 
from its holder, cooled in water, and replaced in a 
inverted position with the thorium iodide deposi 
again at the top of the tube. For many of th 
photographs, this 3-min running time was sufficient. 
for longer exposures it was necessary to cool and 
invert the tube several times. This change could be 


sceCc. 


made within 30 see, 
The primary standard on which the thormum 
measurements were based was a NBS-Megger 


He'* lamp [7], that contained 1 mg of Hg" and 
argon at a pressure of approximately 3 mm of Hg 
The mercury tube was excited by a second dit 
thermy generator at the same frequency of 2,450 Me 
The Hg" source was cooled by an air stream froma 
small blower, 
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3. Experimental Method 


The optical arrangement was the same as that 
ysed by the authors in the previous investigation 
of iron wavelengths from an electrodeless tube [1]. 
For this reason, only a brief description will be given 
here. The Fabry-Perot interferometer consisted of 
aluminized quartz plates separated by 3 sets of 
invar pins either 25, 40, or 50 mm in length. Good 
interference fringes were also obtained with plate 
separation of 67.5 mm, but lack of time did not per- 
mit the exploitation of the higher orders of inter- 
ference. 

The evacuated 10° mm Hg) interferometer 
was mounted externally in parallel light. Both 
thorium and Hg'* lamps were placed on the optical 
axis but at right angles to it and parallel to the 
spectrograph slit, so that both sources were ob- 
served side-on. The Hg'* lamp was farther from 
the slit so that it could be imaged inside the thorium 
tube. From this point on, light from both lamps 
traveled the same path through condensing lens, 
interferometer, projection lens, and thence through 
the slit and spectrograph. The quartz-fluorite 
achromatic projection lens had the proper focal 
length (50 em or 100 em) to form six or more circular 
fringes on the spectrograph slit, which was 25 mm 
long and either 0.15 or 0.30 mm wide. Accuracy 
of measurement was thus favored by the large scale 
of the interference patterns and by a slit width that 
made the fringes approximate rectangles. 

Dispersion in the horizontal direction was _ pro- 
duced by a 15,000-line-per-inch concave grating in 
the Wadsworth mounting. The curvature of the 
plate holder was decreased so as to favor the hori- 
zontal focus and thus produce sharp images of the 
Fabry-Perot fringes from both lamps, simultane- 
ously over as large a range as possible. This in- 
strument was used to photograph diametral sections 
of interference patterns for waves between 4400 and 
7000 A, with a slit width of 0.30 mm and a reciprocal 
dispersion of 5 . mm in the first-order grating spee- 
trum. Because of the higher density of lines in the 
ultraviolet spectrum of thorium, the region from 
$275 to 4500 A was photographed in the second-order 
grating spectrum with a reciprocal dispersion of 
245 A/mm and slit width of 0.15 mm. Under these 
conditions most of the measured interference patterns 
were free from overlaps of neighboring lines, but if 
many more thorium wavelengths are to be measured 
larger grating dispersion or smaller slit width will be 
necessarv. In each wavelength range several spec- 
trograms were made with different exposure times, 
and with slight mechanical-pressure changes in the 
interferometer spacing to alter the interference 
configurations or fractional orders for all spectral 
lines, 

The interferometer plates were the same ones used 
previously for measurements of iron and mercury 
wavelengths [1, 12, 13]. The last series of measure- 
ments [13] suggested that the dispersion of phase 
change in the aluminum films was very small between 
3300 and 6000 A, so no corrections on account of 


| 


| 
| 


lengths. No significant differences were found be 
tween results from the 25-, 40-, and 50-mm etalons. 

All mercury lines that appeared with sufficient 
intensity were measured along with the thorium 
lines. The green (5462.2705 A) and violet (4047.7144 
A) lines of Hg"* were chosen as standards respectively 
for the long- and short-wave interference spectro- 
grams described above. The values of other Hg" 
lines measured on these spectrograms agreed quite 
well with unpublished wavelengths [13], the greatest 
difference being 0.0002 A. This supports our as- 
sumption that any possible corrections to these 
thorium preliminary wavelengths on account of 
dispersion of phase-change in the interferometer 
films are negligible. 

The interferograms were measured at the National 
Bureau of Standards with the interference comparator 
designed by K. Burns [14], and at Purdue University 
with a comparator made by Carl Leiss. In order to 
measure as many thorium wavelengths as possible 
within a given time, it was found desirable to measure 
just 5 diameters of each interference pattern. As 
mentioned before [7], by measuring the diameters of 
5 interference rings, the least-squares evaluation of 
the fractional order at the center may be made with a 
minimum of arithmetical computation because some 
of the steps are easily performed mentally. Even 
with this conservation of time and labor the thorium 
wavelengths presented in this paper required the 
bisection of more than 12,000 interference fringes and 
calculations for more than 1,200 individual patterns. 

Another way to reduce the time and labor involved 
in interferometric comparisons of wavelengths is to 
begin with the largest orders of interference consist- 
ent with finding the correct integral order for each 
line. If grating measurements of wavelength are to 
be refined and replaced by interferometric values, the 
initial orders of interference should not be so great 
that errors in the grating values lead to wrong in- 
tegral orders. For example, if a thorium wave- 
length near 5000 A is to be measured relative to 
5462.2705 A of Hg’, the whole order of the thorium 
wavelength with a 25-mm etalon: will be about 
100,000, and the error in the grating value must be 
less than +-0.02 A if it is to yield the correct integral 
order. We assumed that the available grating values 
of thorium wavelengths justified beginning the inter- 
ferometric observations with a 25-mm etalon and then 
testing these values with a 40-mm etalon. If the 
integral orders of interference are in error because of 
too large errors in the grating values, the 25- and 40- 
mm etalon values will differ radically, but if the inte- 
gral orders are correct, the results will agree within 
a very small possible error of observation. This bold 
procedure saved much time and labor by eliminating 
interferometric observations with all smaller etalons 
which do not contribute anything to the accuracy 
attainable with the larger etalons. 

Two independent sets of grating values of thorium 
wavelengths were available, one in the M/7 Wave- 
length Tables [2] representing wavelengths emitted 
by a thorium are, and the second representing un- 
published wavelengths from thorium halide lamps 


phase change were applied to these thorium wave- ! measured by R. Zalubas in the National Bureau of 
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TABLE I. Thorium wavelenath 
ie Ava 
MIT 1165. 813 1166. 
NBS 1165. 782 1166. 
MIT 3869. 639 3870. 
NBS 3869. 665 3870. 
MIT 3854. 547 3855 
NBS 3854. 512 3855 
MIT 3421. 189 3422. 
NBS 3421. 215 3422. 


Standards. Both sets were referred to the adopted 
secondary standards of wavelength obtained from the 
iron are at atmospheric pressure. We gave preference 
to the MIT grating values because these were ob- 
jectively measured and calculated by automatic ma- 
chines [2], and therefore entirely free from human 
errors. This was a mistake; the only errors in inte- 
gral orders of 25-mm etalons occurred because of 
errors larger than 0.02 A in the MIT grating values 
as shown in table The first line shows in succes- 
sion the MIT wavelength in air, its value in vacuum 
[15], and discordant results from 25-mm etalons and 
from 40-mm etalons. The second line gives similar data 
for NBS grating wavelengths, but the results from 25- 
and 40-mm etalons are in perfect agreement, because 
the integral orders were corrected as indicated -by 
the NBS grating values. The interferometric values 
for these four lines show that the MIT grating val- 
ues have errors of + 0.047, 0.025, +0.036, and 

0.021 A, respectively, whereas the NBS errors are 
-0.016, +0.001, +0.001, and +0.005 A. No ex- 
planation of the large and erratic errors in the MIT 
grating values is offered, but this experience should 
be a warning that some 7-figure MIT wavelengths 
may not be trustworthy beyond 5 figures. 


4. Results 


The thorium wavelengths obtained as a result of 
this investigation are presented in table 2. Since 
the interferometer was evacuated and the separation 
of the interferometer plates was always expressed 
in terms of vacuum wavelengths of Hg" standards, 
all we measured wavelengths of thorium are valid 


only in vacuum. Column | displays only that part 
of bens vacuum W fore asc that can be expressed in 
whole angstroms, the fraction that must be added 
to this is given in columns 2, 3, and 4, depending on 


it was measured respectively with 25-, 40-, 
r50-mm etalons. <A figure in parentheses following 
each fraction shows the number individual 
terminations. Column 5 contains the average pre- 
liminary vacuum wavelength resulting from all the 
observations. In column 6 the vacuum wavelenghts 
of column 5 are changed to values in standard air 
with the aid of Edlén’s conversion table [15]. 

umn quotes unpublished estimated relative 
tensities of thorium lines as given by Zalubas, 


whether 
> 


of de- 


(‘ol- 
in- 
and 


‘ 





grating and interferometer values 


8 


Avae (25 mm) Change Avae (40 mm 
in order 
QRZ 1166. 9750 1166, OS32 
956 1166. 9403 | 1166. 9403 ' 
736 3870. 7308 3870. 7416 
762 3870. 7606 l 3870. 7604 
} 
5. 640 3855. 6333 3855. 6222 ' 
5. 605 3855. 6036 l 3855. 6035 
170 3422. 1673 3422. 1618 
196 3422. 1908 l 3422. 1910 


column 8 indicates the particular 
each wavele ngth be longs. 

It is a common characteristic of all « ‘complex spectra | 
that neutral atoms dominate the visible, and singly. 
ionized atoms predominate in the ultraviolet, jp 
these preliminary measurements of thorium waye. 
lengths with Fabry-Perot etalons no difference jp 
character could be observed between the Th1 and | 
Th 1 lines, and there is no reason to suspect that one | 
category is more accurately measured than the other, | 
Unfortunately none of the thorium wavelengths has 
been measured a sufficient number of times with any 
given etalon to compute a meaningful probable errr 
or standard deviation. A considerable number of 
the wavelengths from a given etalon are limited tog 
single observation and a comparison of these with 
the over-all averages may give an impression of the 
uncertainty in individual observations. 

Because of the lack of sufficient observational 
material to calculate probable errors of the thorium | 
wavelengths, an attempt was made to test their 
relative values by means of the combination principle 
Although the Thi spectrum has been extensively 
analyzed [16] no recurring wavenumber differences} 
were found between Thit lines in table 2. Some} 
what better fortune was encountered with Thi 
lines. Even though the published analysis of Thi 
extremely fragmentary [17] it now being 
extended by Zalubas [18] who has reported that 27 
of our interferometric values of Thi lines test. the 
constancy of 4 wavenumber differences, viz., 2869.26], 
3687.990, 5563.143, and 6362.396 K (=em—') between 
low-energy levels established for neutral thorium 
atoms. These are presented in table 3 as 14 pails 
of vacuum wavelengths in A, and their reciprocals, 
or vacuum wavenumbers in K. The mean deviation 
of the 14 differences from their arithmetical averages 


spectrum to i 


is Is 


is 0.001, K, that is 1 part in about 20 million. — The 
pairs numbered 2 and 3 account for two-thirds 
of the summed deviations; if these are ignored the 
remaining 12 exhibit an average departure of only 
0.0005 K, or 1 part in about 40 million, but this 
may be discounted as accidental. 

The 222 wavelengths (3288.7356 to 6991.5839 A 


of thorium radiations in table 2 are presented 4% 
a preliminary set of improved standard wavelengths 
for the measurement of highly-dispersed complex 
The aver interval between these stand-| 


spectra. age 


98 





PABLE 2 
a. W 
avele 
ngths 
in thori 
um & 
spectre 
1 


~ 
“I 
2 4 





i 
Whol 
le Frac 
act 
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TABLE 2 Wavelengths in thorium s pectra—Continued 


Fractional angstrom 


Whole Vacuum Al Relative Spectrum } 


angstrom Intensity 
25 mm(obs 10 mmi(obs 50 mm(obs 

5462 2705 (5 2705 (5 2705 (1 5462. 2705 D460. 7530 Hy } 
5453 7338 (1 7343 (2 5453. 7341 5452. 2188 250 I 

5452 6215 (4 O2Z11 (2 6216 (1 5432. 6212 5431. 1116 300 I 

5427 1859 (2 IS65 (3 5427. 1863 5425. 6781 250 I] 

5418 9915 (5 9916 (5 9916 (1 5418. 9916 5417. 4856 200 | 

5409 1568 | 1570 (3 1573 l 5409, 1569 5407. 6535 200 ] 

5388 1093 (1 1083 (1 1OS7 (1 5388. LOST 5386. 6109 00 | 

5345 0672 (3 OG7S8 (3 OOS] l 5345. 0676 5343. 5813 500 | 

5328 1573 (3 75 (3 1572 (1 5328. 4574 5326. 9755 LOO | | 
9278 9691 (2 9687 (2 ‘ 5278. 9689 0277. 5002 100 I] 
5259 $245 (3 S245 } 5259. S245 5258. 3609 300 | 

9232 6158 (3 O160 > 9232. 6159 d231. 1596 YOO | 

5178 1028 (3 1023 (3 5178. 4025 5176. 9606 Loo | 

5160 O408 (3 0413 (3 5160. O411 S1SS. GO4I 700 | 

5155 O78] | 6789 (3 9155. G7S87 5154. 2429 Loo | 

5116 1694 (2 1700 (2 5116. 4697 5115. 0443 250 I i 
5069 S865 (2 3870 (2 5069. 3868 5067. O739 G00 | 

5051 2039 (2 2039 (2 D051. 2039 D0AD, 7959 100 I] 

5030 O5S81 (2 0594 (2 5030. OSS88 5028. 6564 100) II 

SOS 6537 4 6542 2 SOLS 6535 SOT 2545 500 II 

5003 1922 (3 1922 (2 5003. 4922 5002. O9OGBS8 10 | 

194] 0205 (4 0204 (2 1941. O2POS 1939, 6417 350 I 

192] IS94 (3 ISS6 (2 1921. IS90O VIO, S155 600 II } 
1896 S214 (4 3216 (2 : 1896. 3215 1894. 9546 350 | 

ISSO O950 (2 0965 (4 ISSO. O96 i878. 733 200 | 
1866 8360 (2 S360 (1 1S66, S360 IS6O5. 4769 350 | 

IS64 5306 (3 5309 (1 tSO4. S307 1S6H5. 1722 L000 I] } 
1842 1949 (2 1953 (2 i842. 1951 1S40. S426 Loo | | 
1809 775 (3 770 (2 1809. 4773 tSO8. 1334 50 | 

1790 7254 (3 7259 (2 1790. 7256 1789. 3867 300 | 

1767 9328 (2 Q334 | 1767. O330 1766. GOO] 200 | 

1753 7427 (2 7432 (2 1753. 7430 1752. 4139 500 I] 

1705 3059 (3 S060 (2 _ 1705. 3060 1703. OS9O7 500 | 

1687 5062 (2 5055 (4 1687. 5060 1686, 1944 1200 | 

1674 9691 (2 G6S9 (2 1674. 9690 1673. 6608 600 | 

1669 TST (3 789 (2 1665. 4788 1668. 1720 700 | 
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1466 5936 (3 DUA0 } 1466, 5038 1165. 3406 on I] 
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1446 Hdb2 2 G1 2 L446. 5561 1445 BOSD tM) | 

1434 2074 ! 2076 (2 $434. POT75 1432. 9628 Hoo [| 

1410 1211 (5 1210 (3 HO. 1211 L408, S828 600 | 

1404 1635 (5 L641 ) 1104. L637 HO?. O70 10 | 

1402 S170 (1 SISS (1 1402. SIS] 1401. 5SS17 1 | 

102 440 (7 440 (2 192. 3440 MOL 1104 Ooo i] 

1383 OOLT (5 OO15 (2 383. OO16 i881. S604 OO) | j 
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TABLE 2 Wavele ngths in thorium spectra Continued 





1 2 3 1 5 6 7 8 
Fractional angstrom 
"um j Whok Vacuum Air Relative Spectrum 
angstrom intensity 
25 mm(obs 10 mm(obs 50 mm (obs 
1379 1069 (5 72 (2 1379. 4070 1378. 1768 500 I 
1375 3535 (2 3538 (4 1375. 3536 1374. 1244 600 I 
1367 1572 (5 1575 (2 1367. 1573 1365. 93038 600 I 
I] 1359 5624 (4 5624 (3 1359. 5624 1358. 3374 H gi% 
1343 1765 (4 1760 (2 1343. 4763 $342. 2555 300 I 
1332 0620 (3 O617 (1 1332. 0619 1330. 844] 300 I 
1310 62302 (4 6307 (1 1319. 6305 1318. 4160 700 ] 
j 1316 1681 (3 1681 (2 1316. 468] 1315. 2544 100 I 
1308 S78 (3 3878 (2 1308. 3878 1507. 1763 700 I 
I] j 130] 0490 (3 O486 (1 101. O4AR9 $299, S393 600 | 
’ 
1993 0176 (3 O178 (2 1293. 0177 1291. 8102 100 I 
1278 5178 (3 5181 (2 1278. 5179 277. 3142 1200 I] 
1274 5599 (3 5603 (2 1274. 5600 273. 3574 1000 I] 
$1258 6946 (2 6941 (2 $258. 6944 $257. 4959 700 | 
1236 6563 (3 6560 (1 1236. 6562 235. 4635 600 l 
123] 6183 (3 6180 (1 $1231. 6182 230. 4268 600 I 
116 O152 (2 OGL (] 1216. O156 $214. R983 200 | 
| P10 0766 (3 O762 (2 1210. O764 $208. 8907 3000 I] 
I $104 IYSL (3 1980 (2 1194. 1980 1193. 0165 900 I 
1179 2375 (3 237 2 1179. 237 L178. 0598 3000 Il 
1166 9403 (3 9403 (2 1166. 9408 1165. 7659 1000 l 
1159 7078 (3 7074 (2 1159. 7076 1158. 5351 S00 ] 
j 1151 1566 (3 1570 (2 H151. 1568 1149. 9865 800 I] 
$133 9191 (3 9191 (1 1133. 919] 1132. 7534 600 I] 
} 1128 5761 (3 5759 (2 $128. 5760 1127. 4117 1400 | 
1116 9199 (3 9201 (2 $116. 9200 $115. 7587 800 I 
} 1109 5788 4 5aT9O] l £109, 5789 L108. 4195 SOO I] 
110} 1985 (3 1984 (2 H1OL. 4984 £100. 3412 1100 I 
10905 9025 (2 9031 (2 1095. 9O28 1094. 7470 1600 II 
1087 6743 (2 6739 (2 1087. 6741 1086. 5205 1600 I] 
1068 49938 (3 5992 (2 1068. 5902 1067. 4507 100 | 
1060 900 (2 3991 (2 1060, 3900 1059, 2525 1000 | 
1047 7144 (3 7144 (2 1047. 7144 1046. 5712 Hgi% 
1044 5366 (3 D371 (2 1044, 5368 1043. 3945 S00 I 
1037 ISSO (2 IST8 (2 1037. 1879 1036. 0475 LSOO a | 
1020 26190 (4 2649 (2 1020. 2649 1019, 1289 1000 II 
1013 6295 (3 62902 (2 1013. 6298 1012. 4950 2000 I 
10090 3433 (3 3438 (2 L009, 3435 1008. 2104 1600 | 
3905 6786 (3 6787 (2 3995. 6786 3994. 5490 1200 I 
398] 2147 (3 2153 (2 981. 2150 3980. 0892 1100 I 
S068 5145 (3 5142 (2 968. 5144 3967. 3919 2000 I 
S950 OS14 (2 OS10 (1 3950. O13 3948, 9636 1000 I 
S934 0243 (3 0243 (1 5934. 0243 3932. 9108 1400 ] 
5924 9103 (3 9106 (1 3924. 9104 3923. 7993 100 I 
S906 2992 (2 2926 (1 3906, 2924 3905. 1861 1500 I] 
ISTO 7606 (2 7604 (2 3870. 7605 3869. 6635 600 | 
3864 5OOS (2 5OLL (2 3864. 5009 3863. 4055 1200 I] 
S855 6036 (2 6035 (2 3855. 6036 3854. 5105 1200 I] 
3843 O497 (2 0499 (2 3843. 0498 3841. 9600 1200 I] 
S40 7829 (3 (837 (2 3840. 7833 3839. 6941 2500 I 
38209 1707 (3 1709 (2 i829. 4708 3828. 3845 3200 | 
381¢ THO] 7692 (1 3819. 7692 3818. 6855 500 I 
S14 1497 (3 1497 (2 3814. 1497 3813. 0674 1200 I] 
3804 1545 (3 1550 (2 SSO4. 1547 3803. 0750 1000 I 
j STR6 6748 (1 6749 (] 3786. 6749 « 3785. 5907 1000 II 
; 
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TABLE 2. Wavelengths in thorium spectra—Continued 
l 2 3 » 6 7 
Fractional angstrom 
Whole Vacuum \ir telative 


angstrom 








itensity 





25 mm(obs 10 mm(obs 
3782 O398 (1 0405 3782. 0402 3780. 9663 350 
3712 HIS (3 1418 3772. 4418 3771. 37038 1500 
3764 OO37 (3 0027 3764. 0032 3762. 9345 1200 
3753 6340 3 6353 3793. 6351 3752. 5685 3500 
3743 9871 (3 9873 3743. 9872 3742. 9231 1100 
3728 9626 (3 9622 3728. 9624 3727. 9022 SOO 
3720 1924 (35 1926 37 20. 4925 3719. 4345 3000 
3712 3595 (2 3597 3712. 3596 3711. 3037 600 
3702 O310 (3 OS13 3702. O83 12 3700. 9780 300 
3693 G171 (3 6171 3693. 6171 3692. 5661 1200 
3683 5343 (3 D347 S683. 5345 3682. 486] L000 
3671 O37 (1 O14] 36071. OL39 3669, GOST 750 
3669 IS42 (5 IS44 3669. 1843 S668, 1396 1000 
3657 72D2 (3 7354 3657. 7353 3656. 6936 1000 
3651 1962 2 1971 S651. L967 3650. L566 
5643 2866 (3 2R67 3643. 2867 3642. 2487 2?00 
5633 8656 (3 S654 3033. S655 3632. S299 1000 
5623 SIR? (3 S20 3625. S281 3622. 7951 SOO 
sO16 1633 | O34 S616. L634 S615. 1324 POOUD 
3613 575 (3 1573 3013. 4574 3612. 4271 1400 
599 1459 (3 1464 3599, 1462 3598, 1196 2000 
3993 SO40 (3 S042 3593. SO41 3592. 7788 2000 
SOS85 1985 (2 19S] S585. 1983 3584. 1753 SOO 
ye 5786 (3 2782 44. DER4 3976. 5573 1000 
S568 PRIS (3 IRIS 3568. PRY 3567. 2635 1200 
S560 1658 (1 1657 S560. 4657 1400 2500 
3552 Hol ) H15OS $552. 4159 555 1O13 LOo0 
545 0306 (3 O300 3545. 0303 jodd. OLF6 1500 
3540 5OSO (3 5OSS S540. 5SOS2 3539. 5SS67 1000 
3519 1094 (3 1095 3519. 4094 SOIS8. 4083 1O0o0 
O12 1610 (3 1613 so12. 1612 Soll. 1570 1000 
3504 7SSO (1 7870 3504. T7875 5503. 7852 500 
3499 O215 4 O216 5499, G2I1G 498, GLZOG Coo 
3494 5174 (2 5174 404. 5174 $493. 5177 2000 I] 
3480 1680 (3 1688 3480. 1683 5479. 1723 R00 I] rds | 
3469 2126 (3 2125 3469, 2125 468. 2193 2000 1 300 
3463 S419 (2 S417 (2 3463. S418 5462. 8500 2000 II UY 
3452 6910 (3 908 (2 3452. 6909 S451. TOL G00) avel 
5443 5H52 4 5651 (2 3445. 5651 6442. 5STSS SOO) f sta 
5454 9827 (3 O831 (2 3434. QS8V29 5435. GOSS 000 vale 
3422 19OS (3 1910 (2 3422. 1909 3421. 2097 2500 | unce 
S413 YSIS (3 9918 (2 3413. 99IS 3413. O127 1800 I hits, 
S406 5346 (3 5348 (2 406. S347 S405. 5575 1400 | rtlo 
3397 7021 } A022 3 307. TO22 5396. 7273 1400 I stril 
54 OOS3 | OOS] 34. OOS 392. OS48 HoO00 " 
Is | 
S386 5033 ; 5O30 (2 S386. 5082 3385. 5312 S00 re on 
338 R300 (3 8305 (2 S381. $308 3380. S595 900 v col 
33" 9438 (2 9440 (2 3375. 9439 3374. 9746 1600 perhay 
8350 5663 (3 5OGS 3359. 5666 3358. GOL 2000 ; 
3352 S10 (3 1813 (2 3352. 191] 3351 92978 2500 Pease 
5358 S206 } S297 (. 535358. S296 de S697 2000 
333 1 $47 (2 143 $331. 4345 3: 1765 1SOO 
3326 O765 (3 O768 (2 3326. O767 3325. 1201 O00 
3325 TOSS (3 TOSS (2 3325. TOSS 824. 7523 2000 
3310 3174 (3 3180 (2 3310. 3176 3309. 3650 S00 
05 ISO4 4 05. ISO4 3304. 2381 3000 
3293 1685 (3 1OS6 3293. 4685 3292, 5202 3000 I] 
S288 7355 (3 1397 (2 3288. 7356 3287. T7885 1600 I] 
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TABLE 3 Test of Thi wavelengths 


Line pait Wavelength Wavenumber Wavenumber 
difference 
| K K 
1494. 5941 22248. 950 2869 260 
= > 2000. 2 
5160. O411 19379. 690 a 
1375. 3536 22855. 296 9860 O77 
2 5003. 4922 19986. O41 —— 
3840. 7833 206036, 356 2860 27 
a 1316. 4681 23167. O89 Ss ae 
3643 2867 27447 744 IRGO 1 
1068. 5993 24578. 483 a 
. 3509. 1462 27784. 367 2869. 261 
» 1013. 6293 24915. 106 en 
3331. 4345 30017. 099 2869. 261 
0 3683. 5345 27147. 838 ios 
" 3310. 3176 30208. 582 2869. 261 
3657. 7353 27339. 321 ee 
i806, 3215 20423. 495 3687 9R9 
s 5975. 3199 16735. 506 = 
“RO 7 Te ¢ IBD — 
- 382! 1708 26118 26 3687. 989 
14540, 2531 22425. 280 
10 3305. 1894 30255. 452 3687. 993 
3764. 0032 26567. 459 et oe 
, ISO) 3215 20423 105 5563 | 13 
= ~_— . “0 eDdedtbe . 
6729. 3157 14860, 352 4 
1° 1258. 6944 23481. 375 5563. 143 
ie 5580. VOT7 17918. 232 er ere 
‘ 1104. 4984 24381. 333 6362. 396 
” 5549. T7170 ISOLS. 937 sate 
‘ 1037 IS70 24769 716 6362 306 
P 5432. 6212 18407. 320 iia 


rds is 10 A from 3300 to 4300 A, 17 A from 4300 to 
300 A, 23 A from 5300 to 6300 A, and 30 A from 
300 to 7000 A. This distribution according to 
savelength corresponds closely to an average spacing 
(standards at intervals of 72 K on a wavenumber 
vale throughout the range 14300 to 30400 K. 
‘ince spectra are always measured in wavelength 
nits, but interpreted in wavenumber units (pro- 
portional to energy) this is suggested as the ideal 
istribution of secondary standards of w avelength. 
lt is emphasized again that the values in table 2 
re only a preliminary list which should be refined 
comparing still higher orders of interference, and 
xrhaps measuring still more wavelengths to de- 
rease the intervals between standards. 


The present list of 222 thorium wavelengths in 
table 2 may serve for immediate direct measurement 
of complex spectra from 3300 to 7000 A, and by 
photographing these standards in a grating’s second- 
order spectrum and doubling the values, they may 
serve also for the measurement of spectra from 6600 
to 14000 A photographed in the grating’s first order. 

Finally, the most important next step is to ex- 
tend the measurement of wavelengths from thorium- 
halide lamps to shorter waves, that is from 3300 to 
2000 A or less. When this is done the iron are at 
atmospheric pressure may be abandoned as a source 
of standard wavelengths (except for rough measure- 
ments) but still retained as an aid in recognizing the 
thorium standards. Because of the highly complex 
character of thorium spectra, and the lack of striking 
groups of lines or useful maps, it will be convenient 
(and perhaps necessary) to juxtapose the simpler 
and more familiar spectrum of the iron are for pur- 
poses of orientation. 


R. W. Stanley acknowledges the support received 
from a grant of the National Science Foundation. 
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An Electrical-Analog Method for Transient Heat-Flow 
Analysis 


A. F. Robertson and Daniel Gross 


An eleetroniec instrument 
fast-time’”’ 


thermal prototype 


type 


is described which permits solution of transient heat-flow 
problems by use of direct analogy to electrical networks. 
simulating the transient involved in times of the order of 107° that of the 

A photoformer type of input-signal generation is used and this, together 


The instrument is of the so-called 


with a variable-frequency master oscillator, permits flexibility of input-signal waveform and 


time seule 
Sore 
discussion 


1. Introduction 


Laboratories in which studies are performed on the 
jro-resistance behavior of building structural com- 
ponents have been frequently asked to predict or 
stimate the fire performance of structures that have 
In some cases data are available 
In most 


ot been tested. 
pon which useful estimates may be based. 


leases, however, the quality of the resulting estimates 





— — 


' 


' 


slargely dependent on the experience and feeling the 
stimator has for the problem. In an effort to assist 
n placing such estimates on a more exact basis, an 
lectronic device for performing the necessary com- 
putations was designed and constructed at the 
Bureau. 

Fire-resistance tests are performed on portions of 
buildings for the purpose of determining the suita- 
bilitv of the construction in withstanding the effects 
f fire. While mechanical behavior may frequently 
limit the usefulness of the structure in this respect, 
it often happens that thermal transmission is the 
ritical factor in determining its” ability to 
fire. The procedure used during such tests 
\}' requires application the furnace-enclosed 
structure Of a time-varying temperature function 
orresponding to the standard fire exposure. The 
procedure also permits heat losses to occur from the 
nexposed portion of the specimen by both radiation 
nd convection, These conditions make it impracti- 
al to use analytical methods for solution of the 
eat-transfer equation, 

Because of this it seemed desirable to use some 
igh-speed approximation methods for computing 
the thermal behavior of structures exposed to fire. 
Both digital and conventional analog computers 
vere considered and used with a certain degree of 
success for approximate solutions of such problems. 
However, it seemed likely that use of the direct 
malogy between thermal and electrical circuits 
might permit greater flexibilitv in the solution of 
roblems and simplify the “coding”. 

The device constructed is somewhat similar to the 
ne developed by Lawson & MeGuire [2]. This 
kes use of the analogy between electrical and 
thermal circuits directly without the need for a large 


resist 


to 


ssembly of electronic mathematieal-operator units 


I Ires in } cket ! te the terature reference t the end « 
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applications of the instrument to the field of fire research 
of errors includes consideration of those resulting from the lumping of circuit 


are described. The 


which, in conventional computers, are assembled in 
appropriate fashion to solve the mathematically 
expressed relationship describing the heat-transfer 
problem. This direct use of the analogy permits 
building a basic signal-generating and measuring 
unit which is much more compact than the usual 
analog computer. The theory, construction prob- 
lems, operation, and limitations of this device are 
briefly described. 


2. Analogies 


In many engineering investigations of the physical 
behavior of mechanical, electrical, fluid dynamic, 
thermal, and optical systems, the complexity of 
analytical or approximate numerical solutions of the 
mathematical formulation of the problem may be- 
come very difficult and excessively time consuming. 
In the past this has often resulted in the use of 
intuitive design procedures. In such an 
approximate prototype of the system involved was 
constructed and experimental studies made of its 
performance. In many instances it was difficult to 
vary the design parameters without extensive re- 
building of the prototype. Because of this, scientists 
and engineers concerned with the design of such 
svstems have, with increasing frequency, made use 
of analogies. Thus an electrical system may be 
analogous to a mechanical, thermal, or other system 
provided there is a likeness between the two systems. 
Usually such a likeness is not one of appearance or 
form but a resemblance of physical behavior. Thus, 
if it can be shown that the mathematical descriptions 
of the behavior of two different systems are similar, 
then it is usually possible to use one system as a 
means for studying the behavior of the other. The 
ease with which electrical circuits can be assembled 
and their behavior studied have made them par- 
ticularly useful as analogs. 

Analogs have long been used by men, maps perhaps 
represent one of the earliest useful forms. Karehhoff 
[3] in 1845 utilized the analogy between current flow 
in a plate and flux distribution between charged 
parallel-line conductors in a uniform medium to 
determine the capacitance between certain configura- 
tions of such conductors. Since then, hundreds of 
publications have appeared on the subject of analog 


cases 
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techniques. Higgins [4] has compiled a very exten- 
sive bibliography of electroanalogic methods in which 


appear references to almost a thousand original 
papers. Soroka’s paper [5] provides one of the better 


introductions to the great versatility of analog 
techniques in the solution of engineering problems. 
His later book [6] provides a detailed exposition of 
many possible analog methods and techniques. 

Beuken [7] presented one of the early descriptions 
of the use of electrical circuits for analysis of transient 
heat-flow phenomena. Paschkis [8] has constructed 
and made extensive use of a similar but larger device. 
Both of these instruments are of the so-called “long- 
time’ type. This implies that the thermal transients 
are studied by means of electrical circuits requiring 
transient periods of several minutes or more duration. 
This permits the use of potentiometric-type instru- 
ments for data-recording purposes. The major draw- 
back to the use of such systems is the requirement 
for use of rather large capacitors. 

Lawson and MeGuire [2] describe a device of the 
“fast-time’’ type. Here the electrical transient 
completed in a few milliseconds and the problem 
repetitively solved, measurements being made by 
means of cathode-ray presentation of the data. This 
latter device was intended primarily for solution of 
problems quite similar to those considered here. In 
their paper they present a rather comprehensive 
exposition of the basis of the electrical-thermal 
transient analogy, together with a discussion of 
means of application to one-dimensional space 
problems and some discussion of the possible errors 
involved in the use of lumped circuit elements. 
Because of the extensive discussions in the literature 
of the various types of analog systems possible it 
would serve little purpose to repeat the descriptions 
here. A brief indication will, however, be made of 
the theoretical basis for use of electrical models for 
solution of transient heat-flow problems. This 
should not be taken to imply that heat-flow problems 
are the only type for which the instrument to be 
described may be used. It consists essentially of a 
voltage versus time signal-generation unit, means for 
application of the generated signal repetitively to an 
electrical network, and provision for measuring the 
transient voltages which are developed within the 
model. Thus, any electrical network that can be 
used as an analogy of another type of svstem may be 
used with this instrument. 

By definition the rate of heat flow 0Q/of in the 
direction of a one-dimensional thermal gradient is 
given as the product of the conductivity &, the area 
A across which heat flow is taking place, and the 
temperature gradient 06/Oz or, 


dQ/dt 


Is 


kAd@/Or, (1) 


the negative sign indicating that flow takes place in 
the direction of a negative temperature gradient. 
The difference in the heat-flow rates into and out 
of a small elemental section along the path of heat 
flow, or the heat-flow gradient, is equal to the product 
of the volumetric heat capacity ps, rate of tempera- 
ture rise 06/Of, and area over which heat is flowing. 
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rye 
Thus 


do o@ 
Or Of 


{ O# 
ae s ° ‘ 
pn Ot ( 
Combining these two equations we obtain, 


k O60 
ps Ou 4 


O8 
or 


This is the one-dimensional form of the waa 
Fourier equation of heat conduction. 

In a similar manner we may develop the equatio; 
for a series-resistance, shunt-capacitance electrical 





transmission line for the particular case in which thi o— 
capacitance is assumed proportional to the Volum) 

of the line2 From the definition of electrical unit] ‘epuT 
and Kirchoff’s laws, 

— 

. aov 
} — ’ ‘ 
r Or ‘ 

_— 
where ¢ is the current, 7 the resistivity of the ling = 
material, a the cross-sectional area and dv/dz q = 
voltage gradient. The current gradient in th — 
direction of flow becomes, 

ey) or 
— C0 ’ 5 
or Or he 
SIGNAL 
; ) INPUT 
where ¢ is the capacitance per unit volume and 
Ov/Of is the rate of voltage rise. Equations (4) and(j} $7 
may be combined to form, 
or 1 Or I 
“y h 
Of reoz 
the t 


Equations (3) and (6) are obviously similar anéj MP 
thus justify the use of the electrical system as ul resist 
analog to the thermal system. In a similar way 

eq (1) and (4) and (2) and (5) may be considered 
analogous. Table 1 has been prepared to list th 
properties of the two systems in analogous form 
assuming equivalence of the cross-sectional areas 
of the thermal and electrical circuits. 

Because of difficulties of constructing within small 
dimensions, electrical models of such transmissiot 
lines in which the capacitance is distributed along 
its length, it is common practice to approximate 
this by a lumped system network. Three sue 
networks, incorporating a series of either T, L, or Pj 
sections, such as those bounded by dotted brackets 
are shown in figure 1. It may be seen that althougl 
the three lines consist of different arrangements @ 
resistance and capacitance elements, the lines an 
identical except for portions of the first and las 
sections. The resistors, R, of these lines thet 
correspond to the total resistance of the thermal 
circuit divided by the number of sections, N, used 
in the electrical model. In a similar way the capact 
tors correspond to the sectional heat capacity 


y, al 
conve 
end 0 
for u 
and } 
Ag 
and 
show! 
eect 
the e 


as a 
will 
Thus 


capac 


simil; 


and f 


? This differs from the geometrical considerations commonly used in determi 
ing capacitance. Such treatment is used, however, to permit development of 
the equations shown and to thus maintain true similitude of the thermal-electriea] 
analogy. ' 
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TABLE 1. 


Thermal 




















, Property Symbol Unit 
un, 
Pemperature “” Cc 
Tine te second 
Heat-flow rate IQ/ot: watt 
a Heat capacity ps joule/em *° C 
' Conductivity k watt/em ° C 
| Leneth Tr em 
, lemperature gradient O/T: C/em 
vell-knoxt Rate of te mperature rise O/ot: C/second 
* equatior 
| etic a2 R/2' R/2 RZ R2 RR R/2 RR 
whic 
ich thi | 3 Re 
he Volum c le c € > 
rical | | SMA = T 
nit wpuT 7 ] l Ry Rp 
= ee | hs SEE — = 
. - HEAT LOSS 
(a) T SECTION LINE TERMINATION 
‘ 
4 R rR ' R , R R 
. 3 Rs 
Mf the lind gemma j¢ le) Le c je j j 
Or/Or INPUT) =—s TT" . eae 5 > 
.t the | Ry, Rp 
1 the —— —— — 
— = HEAT LOSS 
TERMINATION 
(b) L SECTION LINE ’ 
 .  o R “ 
5 : T Me 
SIGNAL  S. _ ,% = ? 
| INPUT S 2 
lume and Fe Fe Fe [°2 |e [er fe felre | Kn, 3, 
4) - L i nt, i> en a oo Oe atte 
and (5 _” HEAT LOSS 
| (c) Pi SECTION LINE TERMINATION 
FIGuRE | Comparison of T, L. and Pi section lines. 


0) 

the thermal model, psz/N. The termination, Ry. 
comprising a varistor, or positive voltage coefficient 
esistor R,, together with series and parallel resistors, 
R, and R,, if necessary, is used to simulate the 
onvective and radiative heat losses from the cool 
ad of the thermal structure. The detailed technique 
for using such elements was described by Lawson 
and MeGuire {2}. 

Again considering the similarity between eq (3) 
and (6) together with the analogous properties 
shown in table 1, it is common practice to substitute 


em as al 


o list the 


nal areas 


hin small 
nsmMissiol 


ted alone dectrical for thermal properties directly and to solve 
ct . . . . 
sroximate the electrical problem, but to think of the solution 
‘ree such 8a thermal one. In performing this conversion it 
' L, orRi Will be found convenient to use scaling factors. | 
brackets Thus rather than substituting directly electrical | 
although “pacitance for heat capacity we have, 
‘ments 0 ps ps 
lines ar c or n , 
and last ” ‘ 
nes thet smilarly for resistivity | 
thermal 
’ sed l ] 
N, ust r= or m=>=;—? 
1e Capac mk kr | 
racity oO , 
pach’y “wd for time 
1 in determin 
velopment ¢ } t, t 
rmal-electrica} t, t ’ 
i mnob mnt, 





Analogous elements of thermal and electrical systems 


Electrical 


Property Symbol Unit 
Voltage t volt 
Time t. second 
Current i microampere 
Capacitance ¢c microfarad/em 3 
Conductivity l/r 1/megohm cm 
Length Te em 
Voltage gradient OP/OT. volt/em 
Rate of voltage rise Ob/ot. volt/second 


where n, m, and @ are the capacity, resistivity, and 
time-sealing factors, respectively. Factors such 
as these may be used both to obtain control over 
the size of the electrical components used as well 
as to provide a legitimate manner for conversion 
from thermal to electrical units as implied by 
these equations. This latter objective may, how- 
ever, also be obtained by the assumption of such 
conversion ratios without their explicit mention. 
Table 2 shows the modifications of the analogous 
thermal and electrical elements when capacity, 
resistivity, and time-scaling factors are employed. 
The relationships developed maintain the identity 
of eq (4) and (5). 


Modification of analogy relationships by application 
of scaling factors 


TABLE 2. 


Electrical system 


Property Thermal 


system Direct analogy Modified analogy 
(m=n=@=1) (m= /kr, n=ps/c, 
o=t,/mnt,) 

Resistivity k r mr 
Capacity Ds h ne 
Time ty t. mnot, 
Length r Vor, 
Temperature 0 t t 
Heat fiow OQV/or: i mv 
Heat-loss resistance Wh R mR; 


As an example of the manner in which an electrical 
model may be constructed it will be assumed that 
we wish to represent a concrete floor slab which is 
10.16 em (4 in.) thick. We will assume that the 
model is to have 20 sections of the L type. This 
problem is shown in figure 2. 

For a resistivity scaling factor of m we have for 
the resistance, 


Lr mre 


kN N 


mR, orR 


r 
Nmk 


Inserting numerical values and using a consistent 
set of dimensions for both electrical and thermal 
systems, 


r 10.16 32.6 
‘"Nmk (20) (0.0156) (m)m 


If m=652, each network resistor is 0.050 meg, a 
convenient value. 
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CONCRETE SLAB 
FIRE 
EXPOSED HEATLOSS 
SURFACE #(@) 
G=t(t) 
~ X=10.16cm( 4in.) - | 
42  =0.508cm(0.20in.) 
L 20 = 
a 
R R R R R 
Se ow wr? 1 
SR 
A L = A i I 7 
ar T= — “>. —y T > 
C C C [c C C 
= L L BR, 
ELECTRICAL MODEL(L SECTION) 
Br URI 7 Vethod ol development ot electrical mode for a ftre- 
exposed concrete slab 
Thermal properties: k=0.0156 w/em ° C (0.90 Btu/hr ft F); p=2.24 g/cm 
140 lb/ft 0.837 joule’s C (02 Btu'lb } Electrica iwlues: R=AwkK 
R,=400 K; R,=varistor; C= 449 pul 


In a similar fashion, for a capacity scaling factor 
of n we have. 


psr ner 0 cy Pst 
= = Hi or -e 
N NA A fi 
Again inserting numerical values, 
esr (2.24)(0.837)(10.16) 0.953 
Nn (PO) (n TD 


If n=2120, each network capacitor is 0.000449 uf. 

Finally, for ¢=1, 

t mut, 652) (2120) Lae x 36° 2... 

Thus, in using the scale ratios proposed, the ther- 
mal time seale has been compressed by a factor of 
about 10°. To provide for a thermal period of 8 hr, 
we represent this by 81080 divisions of 6 min 
each. Each division therefore to 
6 60/1.38< 10°=261 psec, and the master oscillator 
Is adjusted Lo operate at a frequency of 3.83 ke. 
The appropriate varistor assembly is connected at 
the end of the network corresponding to the upper 
surface of the floor slab to represent radiative and 


corresponds 


convective heat losses. 

If, as is often likely, it is found that capacitors or 
resistors of the particular sizes suggested are not 
readily available the value of m and or 7 may be 
changed and the frequency similarly adjusted. The 
ease with which this may be done is one of the de- 
sirable features of this particular instrument It 
not only greatly facilitates construction of models 
but, as shown in table 2, permits the use of a given 
model to simulate a whole range of structuré sizes 





simply by independent change of the time seale ; , 
a change in @ When this is done it is important 4 
remember that all dimensions of the thermal svster 
are changed. Thus in the analog of a multicom 
ponent thermal bod, a change of time Stale wil 
result in the simulation of a new thermal System jy! 
which all dimensions in the direction of heat flow are 
scaled in the same ratio. If special terminationd 
have been used to simulate heat-loss conditions for 
such a system, care must be taken to adjust theg 
separately according to the sealing relationships 
shown in table 2. i. 


3. General Description of Instrument 


The instrument constructed of the “fast. 
time’? type in which a particular problem is solyed 
repetitively at a speed much higher than that used 
for the thermal prototype. The ratio of the speed! 
of the electrical operation to the speed of the ther- 
mal process was of the order of 10° This requires 
cathode-ray presentation for read-out of the repeti- 
tively solved problem. In order to reduce the cost 
and time involved in its construction, commerggl 
electronic instruments were used wherever it seemed 
practical. 

A block diagram of the device is shown in figureg, 
This has been divided into three interrelated fune 
tional units: (a) The signal generator; (b) the elgg} 
trical model; and (c) the measuring circuits. The 
three power supplies are not shown. 

The signal generator unit selected was of the 
photoformer type [9, 10]. This type was chosen 
because of the great flexibility of waveforms possible 
with such a system and the ease with which they 
could interchanged. A cathode-ray tube js} 
arranged with an opaque mask covering a portion 
of its face. A phototube with appropriate amplifier | 
is arranged to view the face of the cathode-ray | 


was 








The amplified outpul of the phototube Is applied to 
the vertical deflection plates of the cathode-ray tube 
in such a way that when the spot is ‘‘seen”’ by the 
phototube it is deflected downwards and forced to 
be partially hidden by the mask. ‘The application 
of a relaxation-ty pe timing sweep to the horizontal 
plates of the cathode-ray tube results in the spot 
following the mask outline on the cathode-ray tube 
screen. Since for given operating conditions the 
position of the spot on the screen bears a fixed rela- 
tionship to the voltage applied to the deflection 
plates, a voltage-time signal may be taken from 
them which is controlled by the shape of the mask. 

In practice, a commercial oscilloscope was used In 
conjunction with a photomultiplier tube as_ the 
generator. A special impedance converter was 
built to modify the balanced high-impedance signal 
applied to the deflection plates into a low-impedanee 
single-ended output signal suitable for application 
to the model 

The electrical model was custom-made for each 
Figure 4 shows 
It consists} 


| 
| 


thermal situation being simulated. 
i photograph of one of the models used 
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‘ich they } a series of resistors representing the thermal 


tube is | Tesistance of the prototype and a group of capacitors 
portion | hunting these resistors to ground and representing 
amplifier | ‘he heat capacity of the thermal system. lhe 
“ay tube | S8embled model assumes the form of a_ plug-in 
plied to Mit, provisions being made for twenty discharge 
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points. 
~The measuring unit comprised a number of cir- 
lits which Tah \ be briefly deseribed as follows: 

1. The master oscillator is a variable-frequency 
quare-wave generator used as a_ time-reference 
source, 

2. The frequeney-divider circuit the 
leeade-time markers and submultiples of these for 
snchronizing signals 

3. The temperature and time-reference unit is a 
mixing device which assembles the calibrated time 
vale with a variable d-c voltage for making “‘tem- 
lerature’’ measurements 


generates 


- used in mn oF , 
as the 2. Phe discharge unit is a special circuit used for 
er was ischarging the capacitors used in the model. This 


ischarge process is accomplished every alternate 
weep of the viewing oscilloscope during the period 
f presentation of the time and temperature refer- 
fice signal 

5. The cathode follower is a special high input- 
mpedance cireuit which permits making voltage 
heasurements on the model without 
incal loading of it 
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instrument, 


6. The electronic switch provides a means for 
alternate presentation of the signal from the model 
and that from the time and temperature reference 
unit. 

7. The read-out oscilloscope provides a means for 
study of the electrical behavior of the model and 
permits direct measurement of voltage correspond- 
ing to temperature changes within the model. This 
oscilloscope is used to generate the time-sweep signal 
for the signal generator. Trigger operation is used 
and the sweep signal applied to the signal generator 
is clamped on alternate cycles in such a way that a 
signal is generated only when required for application 
to the model. 

Figure 5 shows the assembled device. The power 
supplies have been mounted below the desk while the 
remainder of the equipment is assembled in two relay 
racks. The electrical model is seen near the mid- 
height of the right-hand rack. 


4. Special Electronic Circuits 


Figure 6 presents the complete circuit diagram for 
the specially constructed portions of the unit. Also 
shown diagrammatically are the commercially avail- 
able electronic units which formed a part of the com- 
plete device. The high- and low-voltage power sup- 
plies were of the regulated type and are not indicated 
in this diagram. The special considerations which 
were involved in the selection of the instruments and 
circuits follow. 


4.1. Signal-Generating Oscilloscope 


The particular oscilloscope used was of the low- 
accelerating voltage, single-anode, electrostatic- 
deflecting tvpe. It was chosen to permit achieve- 
ment of high linearity between mask size and output 
signal. It was fitted with a tube having a short 
persistence type P-15 phosphor. The short phos- 
phorescent decay time, down 90 percent in 2.8 ysec, 
of this phosphor was found necessary to permit 


close following. of transient mask shapes. The 

















photomultiplier tube associated with the oscilloseop, 


was selected to have a spectral Sensitivity which 
permitted efficient use with the cathode-ray tube f 
phosphor. The high-output sensitivity of this tube 

was also a factor involved in its selection. A dark. [— 
ened hood was of course used to shield the photo-| & - 
multiplier from room illumination. In using thjc! 
photoformer circuit it was found that the delay time 

of the oscilloscope direct-coupled amplifiers caused 
considerable hunting. This had the effect of super} — 
imposing a small a-c signal in excess of 80 ke on the Pon 


signal being generated. The presence of this high-| 
. . Boe 2 - . © 
frequency signal was not considered objectional for! peur’ 
normal use in solving problems since it was com.t 











pletely filtered out by the first section of the network 





to which it was applied. It could be greatly reduced 2¢™Me 
in amplitude by shunting the photomultiplier output} !™ 
by a 150 ywuf capacitor. A more proper way off Tlie‘ 
eliminating this hunting would involve design ‘of qj blank 
special amplifier having short time delay and possibly | mast 
incorporation of negative feedback. In_ using aj® "S 
signal generator of this type it is essential that the} ou" 
spot remain on the tube screen during periods when! ™ | 
signals are not being generated. To permit this 















tvpe of operation, the generating oscilloscope was 








modified by insertion of clamp circuits on the hori- | Th 
zontal deflection circuit to prevent the spot from but 0 
leaving the right-hand edge of the screen. The feet 
circuit used is shown in figure 7. It does not ingle 
abruptly limit spot deflection at one position but has} j»ee 
the effect of greatly reducing the deflection sensitivity | ji -e¢ 
y . 9 — when the spot approaches the edge of the screen.} nore 
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tional for FIGURE 7 Clamp for horizontal deflection circuil of generating 


WAS com-| 


» Network! a 7 , a 
v reduced semed desirable to assist in preventing “burning 


er output of the phosphor of the cathode-ray tube screen. 
Way off The oscilloscope used was provided with a beam- 
sign of qy Hanking cireuit for removal of the return trace 
| possibly} image. This feature was undesirable in the unit 
using | used for signal generation and the instrument 
that the ontrols were accordingly set to remove this circuit 
ods when! from operation. 
rmit this 
rOpe was 
the hori. This circuit was designed to combine the balanced 
pot from but out-of-phase voltages applied to the vertical de- 
en. The} feeting plates of the generating oscilloscope into a 
does not ingle-phase output signal with ground as a refer- 
n but has) Capacity coupling was used, although 


Ost illoscope ° 


4.2. Impedance Converter 


nce level. 


NSIUIVILY | fireet coupling if practical would have provided for 
© screel.| nore rapid and easy adjustment of the output level 
sily and) with respect to ground. The pentode is used here 


as a unity-gain phase inverter and the twin triode 
as an averaging circuit. This was acceptable as a 
‘substitute for a true addition circuit because of the 
very similar size of the two signals being averaged. 
The output impedance of the converter as measured 
by the shunt resistance required for 50-percent at- 
tenuation of signal was found to be about 250 
ohms. 





4.3. Sweep Clamp 


The horizontal timing sweep used for forming the 
mput signal was obtained from an internal connec- 
tion to the sweep-generating circuit of the viewing 
wcilloscope. The sweep clamp provided a means 
for rejecting every alternate sweep cycle and thus 
oly generating an input signal when it was necessary 
to apply one to the model. 





4.4. Cathode Follower 
This cireuit is used as a high-impedance probe 
eeding a wide-band impedance converter. The 
arcuit while probably not original with us has the 








INTERN 
| |@ wique feature of permitting high-impedance opera- 
4) tion with large input signals. The potential divider 
3 etween plate and cathode of the first stage of the 
win triode serves as a variable plate-supply source 
__| | br the GAK5 used as cathode-follower input stage. 


ln this manner large voltage signals can be followed 
vhile at the same time limiting the potential drop 
«ross the input tube. The input impedance as 








measured by the decay of voltage on a capacitor was 
found to be at least 8,000 meg. Frequency response 
was highly linear up to 25 ke. The output of this 
circuit is mixed with the time ané calibration signal 
by the electronic switch. 


4.5. Master Oscillator 


The time-reference standard for the unit comprises 
a variable-frequency square-wave oscillator. The 
output of this unit was fed to both the temperature 
and time-reference unit, and the frequency divider. 
The former comprised differentiating, rectifying, 
and mixing circuits. 


4.6. Frequency Divider 


This circuit made use of ‘‘Dekatron”’ glow-counter 
tubes of both scale of 10 and seale of 12 types in 
series. The former permitted accentuation of each 
tenth timing mark while the latter permitted genera- 
tion of a synchronizing signal as a submultiple of 
the scale of 10 count. The output from each of 
the 12 cathodes of this tube were separately avail- 
able so that considerable flexibility was possible in 
generating a time signal appropriate to the input 
signal being used. The output from this seale of 
12 count was used to synchronize the time sweep of 
the viewing oscilloscope and the switching frequency 
of the electronic switch. An internal connection 
was made to the electronic switch to provide a 
square-wave output in synchronism with the trigger 
input signal. This output was used to operate both 
the sweep clamp and the discharge circuit. 


4.7. Discharge Circuit 


This unit comprises three circuits: (a) A twin 
triode amplifier and phase inverter, the output of 
which is split into two signals to drive; (b) the dis- 
charge triodes for each individual section of the 
model; and (c) the clamp-diode reference-level 
generator. This last circuit was necessary to per- 
mit balancing out the potential drop aeross the clamp 
diodes and thus to permit discharge of the network 
model to the ground-reference va In the original 
unit, 20 discharge triodes consisting of 10 twin 
triodes, were used permitting complete discharge of 
a 20-section electrical model. The discharge cir- 
cuit was assembled as a unit in a relay-rack section, 3 
in. high. <A total of 3 discharge circuit assemblies 
could be installed in the relay rack illustrated in 
figure 5. This corresponds to a model requiring 60 
discharge points. 


4.8. Read-Out or Viewing Oscilloscope 


The same type of oscilloscope was used for obser- 
vation of behavior of the model as was used for the 
signal generator. This was done largely because of 
the facilities it provided for driven-sweep operation 
from an external trigger signal. This together with 
the requirement of a highly linear time-sweep source 
for the signal generator are the only essential require- 
ments of the unit. High linearity of the vertical 
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deflection amplifiers is not essential because of the 


method used in) comparing two different input 
signals. , 
Figure 8 presents the waveforms observed with 


the use of an oscilloscope with 2-meg input imped- 
The points on the circuit at which measure- 
ments were made with respect to ground are indi- 
cated in figure 6 by small encircled numbers. These 
waveforms have been presented for the case in which 
the standard time temperature curve was 
being used as the input signal and the master 
oscillator was set for 24-ke operation. 


ance, 


versus 


5. Performance and Limitations 


The performance of the device and the limitations 
on its use will be described in terms of the errors 
involved and a comparison between the results 
obtained by means of analog measurements and 
those obtained analytically and experimentally. To 
date, our main interest: has with one-dimen- 
sional problems, but the extension to 2 and even 3 
dimensions presents practical rather than theoretical 
difficulties. 


been 


5.1. Errors 
The errors involved in the solution of transient 
heat-flow problems by means of the electrical 


analogy described may be divided into (1) errors 
due to “lumping,” also known as truncation errors, 
and (2) errors due to experimental sources 
an iterative process is not involved, the errors due to 
rounding off introduced in numerical computations 
are absent. 


Since 
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o. Errors Due to Lumping 


These are errors introduced when a line CONSistin 
of distributed parameters is represented by one com 
posed of a finite number of sections or lumps, Th 
error is a function of the boundary and initia] 
ditions, the number of lumps, the type of lumpin, 
element used (i. e.. whether . L, i). and thi 
ratio of the time at which measurements are taken ¢ 
the resistance-capacitance product of the electriegt 
model (the Fourier number). Error analyses are g 
far available for only a very few types of problems 
Lawson and McGuire [2] have shown that for a semi 
infinite slab, the surface temperature of which wad 
instantaneously raised to and maintained af 4 
constant level, the temperatures for T, Ly, and p 
section networks are identical after the first sections 
Klein, Touloukian, and Eaton [11] have analyzed 
the distribution of errors in a finite bar (repeeseall 
by different numbers of T-sections) one end of whieh 
is instantaneously raised to a temperature that js 
maintained constant thereafter while the other end 
is maintained at a constant reference temperature 
They found that departure of the lumped system from 
the distributed system at the center of the bar was 
reduced (after high positive departures during thd 
very early part of the transient) from a maximun 
of —6 percent with a 5-section model to a maximum 
of —2 percent with an 11-section model. Greate 
departures were found toward the leading seetion 
in a lumped system so that this point may vied 
unreliable information for the earliest part of the 
transient. Clarke [12] has analyzed the error inal 
finite bar represented by different numbers of L 
sections, one end of which is supplied with a constant 
heat flux, the other end being perfectly insulated} 
He found that departure of the lumped system from 
the distributed system at the insulated end of the 
bar was always positive and decreased with the time} 
at which measurements are made. Friedmann {ij 
has estimated an upper bound for the global or max-| 
mum error between the exact solution of the heat 
conduction equation examined by Klein et al. and 
the exact solution of its lumped parameter analog 

The numerical solution for the distribution of 
lumping errors in the problem of a finite bar with 
perfectly insulated sides, one end of which ts instan- 
taneously raised to a temperature that is maintained 
constant thereafter while the other end remains per- 
fecthy insulated, seems to have been overlooked by 
workers in this field. The solution for the lumped 
system of T-sections was found in Jaeger [14], while 
the exact thermal solution was given by Carslaw and 
[15]. Some of the results of the analytical 
calculations are shown in a series of curves in figures 
Qand 10. The comparison between the errors at the 
center in each of the lumped circuits as shown it 
figure 9 indicates the improvement in the results as 
more sections are used to represent the distributed 
svstem For a given lumped circuit of 10 sections 
error curves for three points corresponding to the 
lengths 0.05, 0.25, and 0.75) are shown in figure 10 
The errors are greatest at the first section and atte 
uate rapidly for points further removed from th 
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A similar series of curves were derived from b. Errors Due to Experimental Sources 


source. — . ° 
the data of Klein et al. for an 11-section system and 


gre shown in figure 11. The similarity in the error 
ne distribution for the cases of short- and open-circuit 
ps. Th terminations suggests that changes in the circuit ter- 


These include (a) consistent errors due to the 
measuring-circuit impedance and response, or due to 
leakage and stray currents arising from the network 
; capacitors or the mounting system; and (b) random 
mination do not result in appreciable changes in the | ¢prors in measuring capacitance, resistance, voltage, 
distribution of errors within a lumped line. time, ete. Paschkis and Heisler [16] have shown that 
for the case where leakage is significant, the error 
due to network capacitor or mounting system leakage 
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electrieg : 
sesareal TT TT) decreases as the number of sections of the network 
problems decreases. ‘The versatility of the time scale employed 
ora semil i TN TT | 7 here permits operation with network resistors below 
hich wad 3 LT 1 \ the megohm range and capacitors below the micro- 
ed at 4 $4. ——+s —~ 1} ae | |} farad range for all materials under study, thereby 
4, and PB é Nits iF minimizing stray currents due to the network ele- 
t sections ° \ ments. The parallel impedance introduced by the 
analyzed } discharge-circuit elements and the measuring circuit 
heme | | is of the order of several thousand megohms. The 
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perature FOURIER NUMBER. Nog able for measurements. The over-all sweep and 
stem from, FioeRre 9. Lumping error at the center of a system represented signal-generation linearity of the device as measured 
» bar wag hy 3, 7, and 21 T sections with a linear-rise input function was found to be 
uring thd open-circuit terminatior within the sweep-linearity specification of the oscillo- 
naximun scope alone, namely, “any 10-percent increment is 
naximum within 10 percent of any other 10-percent incre- 
Greater ment ...’. It is estimated that for the one- 
gy section! [| | dimensional case errors due to all experimental 


sources (which limit the accuracy of the method) are 
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finite bar with perfectly insulated sides, one end of 
which is instantaneously raised to a temperature that 
is Maintained constant thereafter while the other end 
remains perfectly insulated. From the preceding 
discussion on errors and since we are using a 20- 
section network representation, it is evident that the 
difference between the analog and analytic solutions 
are principally due to experimental sources. The 
maximum error is seen to be within the estimated 2 
percent of the applied signal. 


5.3. Comparison With Experimental Results 
Figure 13 is a comparison between the analog and 


experimental results on a horizontal 4-in. concrete 
slab the lower surface of which was subjected to the 
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Figure 13 Temperature rise al center-plane of a 4-in 


concrete slab subjecte d to standard fire er posure 


, Analog solution s=0.837 joule/g ° C 2.09 joules/g 


Cc experimental results 


inalog solution s 


standard fire exposure. This problem was set up as 
an illustrative example in section 2 and figure 14 
shows the form of read-out oscilloscope display for 
such a problem. In normal use only two curves 
would be displayed simultaneously; the calibration 
scale, and a temperature versus time curve for the 
position under study. However, for purposes of 
illustrating a larger portion of the solution, time- 
temperature curves for a number of locations within 
the specimen as well as the two surfaces have been 
recorded by means of multiple photographic ex- 
posures. 

Convective and radiative heat loss on the upper 
surface has been represented by the appropriate 
varistor assembly. The analog result was obtained 
using constant values of density, specific heat, and 
thermal conductivity as noted in figure 2. The de- 
vice is not set up to handle temperature or time- 
varying properties although these variations over 
the temperature range considered may be appreci- 
able. Since the values themselves are subject to 












Vultrple -er posure photograph of read-out oscillo. 


Figure 14 
scope display for transient one-dimensional heat flow in @ 


homogeneous slab. 


Large time-scale markers correspond to 1-hr time periods 


considerable uncertainty, their variation with tem-| 
perature may be considered only a_ refinement, 
However, the effect of variable properties may be 
quickly investigated by duplicate tests at several 
different but constant values of the property in ques- 
tion. Or, alternately, by designing a line such that) 
the thermal constants represented at each section 
vary as the mean temperature likely to be obtained 
at those points. The broken line in figure 13 repre 
sents the temperature curve for a specfic heat of 
2.09 joules/g ° C (0.5 Btu/lb ° F) and illustrates how} 
the device may be used to estimate thermal properties 
when direct thermal property measurement is im- 
practical. The extent of the physical changes taking| 
place during fire exposure and the role that moisture; 
plays are a few of the uncertain variables which may | 
cause considerable difference between experimental | 
and analog results. 


6. Conclusions 


An electronic instrument has been constructed 
which permits the solution of transient heat-flow 
problems by means of analogous resistance-capac- 
tance electrical networks. On the basis of the work 
completed, the reliability of the instrument as wel 
as its flexibility with regard to input-signal wave- 
form and time scale have been demonstrated. For 
an assumed one-dimensional problem in which ther- 
mal properties were independent of temperature and 
high accuracy was achieved in the construction of the 
electrical model, the errors were found to be within 
2 percent of the applied signal. The instrument has 
been applied to a variety of problems in the field of 
fire research where a rapid estimation of the thermal 
endurance of structures of given thermal properties 
may be obtained. Although not vet attempted, the 
method may be extended to two- and three-dimer-} 
sional transient heat-flow problems. 
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Instrument and Components 
List 


8. Appendix. 


The special components and instruments used in 
the assembly are listed below. It is quite possible 
that other equally satisfactory components could 
have been used. 


Generating oscilloscope. DuMont type 340 R, cathode ray 
tube type 5AQP-15, short-persistenee phosphor. 

Viewing oscilloscope. DuMont type 340—-R, cathode ray 
tube type 5AQP-1. 

Oscillator. Rutherford Electronics Co., model G1 pulse and 
square wave oscillator. Frequency range 1 to 100,000 pulses 
per second 

Electronic switch. DuMont type 330, arranged for triggered 
operation. 

Multiplier 
response. 

Photomultiplier power supply. ¢ 
1,000 v “ Uniplug”’. 

Counter tubes. Distributed by Atomie Instrument Co. One 
type GC10D single decade (scale of 10) and one type GS12D 
(seale of 12) **Dekatron’’ cold cathode glow transfer counter 
tubes. 

Direct-current voltage supplies. Low-voltage regulated sup- 
ply,-150v, 125 ma. High-voltage regulated supply, + 300 v, 
300 ma. 


phototube. RCA—6328, 9-stage type with S—4 


‘, J. Applegate model 123A, 


WasHInGTon, April 29, 1958. 
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and of spin-orbit coupling on an 
Only fourth-degree terms in the crystal- 
and eigenvectors of the resulting 
704 electronic computer at 
where D is a erystal-field param- 
Perturbations by an external 


magnetic field are computed using second-order perturbation theory. All eigenvalues are 
tabulated for D/i 10-4, 2X 10-4, 10-*, —10°°, —2X10"*, —5X10-3, —107%, 
10 They can be adapted directly to : finite nonzero positive value of ¢ and, by 


interpolation, for values of D/¢ 


calculated and displayed as 


5 102 em 
ot i * 1s discussed 
| 1. Introduction 


The magnetic and optical properties of ionic and 
atomic species in solids are affected considerably by 
| their neighbors in the lattice. Crystal-fiel i theory [ij 
furnishes a simple but reliable technique for taking 
uch effects into account. In this paper are dis- 
ussed the combined effects of the erystal field and 
the ionic spin- orbit —s on the magnetic sus- 
eptibility of ions with an # electron configuration. 
Examples of such systems are the Nd** (4f) and U 
5f) salts: in the discussion the results are applied 
to some of the latter. Penney and Schlapp [2] have 
— the properties of ions with / configuration 
a crystal field of cubic symmetry, with the assump- 
tion that the energy of spin-orbit coupling is very 
arge (effectively infinite) compared to the crystal- 
liield effect. They considered terms of fourth degree 
n the potential expansion, but neglected terms of 
sixth degree. Recently Hutchison and Candela [3] 
have treated the problem taking the sixth-degree 
terms in the potential into account, but they also 
assumed the effect of the ionic spin-orbit coupling 
to be large compared to the erystal-field splittings. 
The calculation reported in this paper makes no 
special assumption about the relative size of the 
erystal-field effect and the spin-orbit coupling, but 
neglects the sixth-degree terms in the potential 
expansion, 
The Hamiltonian of an ion embedded in a crystal 
lattice, in the presence of an external magnetic field, 
may be written schematically in the form 


H H ¢ D-S+-H wl T. (1) 


Here H, is the sum of three tvpes of terms: The 


KInetie energy of the electrons of the ion, the potent ial 


_—— 
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, e United States Government 
Figures in brackets ind te the lit ire reference t the end of this paper 


not tabulated. 
a function of temperature and D/¢, for ¢ 
The relation of the results of the calculations to measurements on compounds 


The magnetie susceptibility of a powder is 
10° em~'!, and 


energy of the electrons in the field of the nucleus, 
and the energy of mutual repulsion of the electrons. 
The term proportional to ¢ gives (approximately) the 
spin-orbit interaction of the electrons. H, sum- 
marizes the perturbations of the electrons arising 
from the neighbors of the ion in the lattice, and p-H 
is the energy of the ~ with magnetic moment yu in 
an external magnetic field 7. The eigenvalues of H 
are required for the calculation of the magnetic sus- 
ceptibility. 


2. Eigenvalues of the Hamiltonian 


The form of the erystal-field energy H, depends 
on the symmetry of the site of the ion. If this sym- 
metry is cubic, H, can be written in the form [1] 
H.=C+- AD, (rt+- yi+ 2{— (3/5) r)+ higher terms. (2) 
Here 


this calculation, x;, y;, 


(is a constant which can be disregarded for 
and 2, are the coordinates of 


i 


the 7th er (the nuc ‘leus of the ion is at the 
origin), 7;= (a?+-y3+- 22), and A gives the magnitude 
and sign of the classical sbi — Two 


cases of cubic symmetry are of interest: The ion 
is at the center of a regular octahedron, with nega- 
tive charges located at the vertices. The ion is 
at the center of a cube, with negative ¢ ae located 
at the vertices. In this calculation only the octa- 
hedral considered, where A= (35ee)/(47r), 
with e the electronic charge, « the charge at a given 
vertex, and r, the distance from center to vertex. 
The ground state of a free ion with f electron 
configuration is *//,, but since the effect of spin-orbit 
coupling is to be taken into account explicitly, the 
erystal-field calculation must involve all 33 states of 
the */7 Russell-Saunders level: *//,, *Hs5, and *7,. 
The —_ convenient representation is SLMsM,, 
with S=1 and L=5. The matrix for H,+¢L-S is 
written out explicitly in this representation, using 
matrix elements calculated from the formulas of 


case 1S 
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Schlapp and Penney [4]. This matrix depends 
on two unknown parameters: ¢, which gives the 
magnitude of the spin-orbit coupling, and /), the 
magnitude of the crystal-field perturbation. J) is 
proportional to A of eq (2), and to (7*), the average 
of r* for the f electrons. These parameters are dis- 


cussed below. The secular determinant of the 
3333 matrix so obtained can be factored into one 
99 and three 88 determinants by rearranging 


rows and columns. The eigenvalues of the matrix 
are computed by solving for the roots of the factored 
secular determinants; the eigenvectors are also ob- 
tained. This computation was carried out on a 704 
electronic computing machine, with a method due 
to Jacobi, adapted for machine use by Givens [5]. 
The eigenvalues are shown in figure | as a function 
of Dit. The eigenvectors (which are also functions 
of D/t) are not shown, because the effect of their 
variation with D/¢ can be taken into account more 
efficiently in ways discussed below. 

The matrix u-// is calculated in the SLA/..M, rep- 
resentation, and then transformed for each value of 


D/¢ to that representation in which H,+¢L-S is 
diagonal, using the eigenvector matrix of H,+¢L-S 


Then the eigen- 


for the appropriate value of Die. 
with second-order 


values of w/7 are calculated 
perturbation theory. 


The eigenvalues of H,+¢L-S+-y-// are of the form 


W=W.+WH+W.P. Here W, is the eigenvalue 
200 ——————————— —_ . —_ 
100 
90 
80 

70 














of H.+¢L-S, and W, and W, give the first-orde 
second-order contributions of uf. 


ho r ang 
_— ‘ . lhe Values ( 
W, f, W’,/8, and W.¢/6" are tabulated in table | f 
the eight values of ))/¢ for which the calculation hal 


been carried out (8 is the Bohr magneton), Eact 
state is designated by a set N of two quantum nu 
bers 7 and 7; 7 is a label for the irreducible represen 
tation of the cubie group [6] of the level W,, and | 
indicates the factor of the secular determinant ly 
which the eigenvalue occurs. The original 339: 
secular determinant factors into four determinants 
(One basis vector is given, to identify the detern; 
nants) 2,: \M,=3, Ms 1), Dy:|—5, —1), Dy:|~4 

1), and Dy:4, 1 D, is a 9X9 determinant, th 
others areS*8. It turns out that the different stated 
belonging to a degenerate level W arise always fron 
different determinants, and further that the set 0 
determinants from which they arise characterizes 
the irreducible representation of the octahedral grouy 
to which the level belongs [6]. Thus nondegenerat, 
levels with } ] belong to 7F,, and those with j 9 
belong to T.. Doubly degenerate levels, slwave ai 
j=1 and 2 belong, of course, to Ts. sty te 





erate levels with j=1, 3, and 4 belong to Ty, those 
with j=2, 3, and 4 to Ts. The method of assign- 
ment is simple. One I, state must be identified as} 
such; the other assignments can be made by inspee- 
tion, by noting whether any nonvanishing matrix 
elements of u-77 connect T, with the unidentified} 
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Calculations are made only for the values shown on abscissa. 
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TABLE 1. Coefficients in expansion of encrgy levels W=W,.4+ W,H-+ WIT? 
Die 10-4 Dit 2xlo-* 
Wee Wie N38 \ Welt W/e8 Wes 
l 6. 0138 0 133.0 oe 6, 0276 0 67.4 
I YS20 0 125.0 a 3 5. 9646 0 61.3 
; YS20) 0). 404 12.4 4,3 5. 646 0.411 6, 22 
t YR20) 404 12.4 4.4 5. 9646 411 —. 22 
! ». YHY2 0 258. 0 3, 1 5. 9104 0 129.0 
2 Y5V2 0 37. 0 3,2 5. 9104 0 18.8 
2 SUuOT 0 37. 31 5.2 ». TROL 0 Ik. 7 
; SO] 1. YOS 12.2 5,3 ». TROL 1. 909 6. 05 
! ». SOOT 1. 9OS 12.2 5,4 ». T801 1. 999 6.05 
1 0. 9936 0 25. 5 4, 1 0. OS72 0 12.4 
; 9936 $. O13 253. 0 4,3 YRT2 3. 006 132. 62 
‘ 9936 $. O18 253. 0 4,4 9872 3. 006 132. 62 
4 9724 0 75. 02 5, 2 0454 0 37.0 
3 9724 2. S87 221.0 5,3 0454 2. 581 117.0 
i Y724 2. S87 221.0 5, 4 W454 2. 581 117.0 
1 ed 0 1400. 0 3, 1 S303 0 789.0 
2 9156 0 75. O06 3,2 S305 0 37.0 
GOSS 0 1520.0 if S107 0 801.0 
7 WO55 2. 506 3L.8 13 R107 2. KOS 15.6 
j GO55 2 wi 31.8 4,4 S107 2. 505 15.6 
l 1. 9024 0 25. 3 3,1 4. UR48 0 —1.34 
2 4. 9024 0 759. 6 3,2 4. 9848 0 2925. 0 
” 4 GORS 0 “41.9 2 4. 9971 0 2509. 0 
; +. VUS5 0. 606 46.4 5,3 4. 9971 0. 688 23. 2 
t 4. 9US5 HO6 6.4 5,4 4. 9971 68S 23.2 
2 5. O168 0 H23. 82 2,2 5. 0336 0 309. 0 
, 5. 0823 0 03, 97 5,2 5. 1652 0 46.7 
' 5. ORS 2. 224 956.0 5, 3 5. 1652 2. 230 178.0 
i 5. OS23 2. 224 056. 0 5 4 5. 1652 2. 230 —478.0 
l 5. 0904 0 1165. 0 4,1 5. 1992 0 635. 0 
2 5 (oo4 0. S86 1000.0 13 5. 1992 0. 590 SOLO 
{ 5. O04 5st} 1000. 0 4,4 5. 1992 590 SOLO 
! 114 0 1190.0 1,1 5. 2291 0 636. 0 
Dit 510-4 D/t 10 
Wee NW/B Wot / se V Wt W/s Wot /e? 
l 6. 0692 0 27.43 i 2 6. 1394 0 14. OS 
1 5. O158 0 25. 66 4,1 5. 8451 0 —13.09 
3 ». 9158 0. 4251 2. 586 1,3 5. 8451 0. 449 —1.37 
4 5. 9158 4251 2. 586 1,4 5. 8451 449 -1.37 
l 5. 8053 0 52. 91 3,1 5. 6317 0 26. 99 
2 5, 8053 0 7. 670 3, 2 5. 6317 0 3. 967 
4 17 0 7.475 b, 2 4. 9024 0 3. 757 
4 107 2. 000 2. 405 5,3 4. 9024 1. 998 1.18 
‘ 4507 2. OOO 2. 405 5, 4 4. 9024 — 1. 998 1.18 
1 0. G6R2 0 4. 692 4,1 0). 9363 0 ~2. 122 
YHS2 2. 978 56.13 41,3 9363 2. 935 — 30. 90 
i gAS2 2. 978 56.13 1,4 9363 2. 935 30. 90 
2 SONS 0 14.09 5,2 7527 0 6. 546 
SHS 2. 578 5. 02 5,3 7527 2. 576 27. 90 
4 SHRS 2. 578 MO. 02 5,4 7527 2. 576 27. 90 
l WS 0 335. 9 3, 1 1203 0 — 198. 0 
2 Wis ) 14.14 3, 2 1203 0 6. 605 
l 5243 0 340.7 4, 1 0421 0 200. 1 
3 5243 2. SOT 6. 128 1,3 0421 2. 492 2.99 
‘ 9243 » I 6. 128 4.4 0421 2. 492 2.99 
l 4. 9628 0 1. 931 3, 1 4. 9280 0 —2. 382 
2 4. O28 0 1137.0 3, 2 4. Y280 0 ~ 561. 5 
2 4. 9936 0 990.3 5,2 4. 9900 0 485.8 
3 4. 9936 0. 6639 9. 030 5,3 4. 9900 0. 621 —4. 33 
; 4. 9936 6639 9. 030 5, 4 4. 9900 621 —4. 33 
2 5. 0840 0 129. 1 = 5. 1680 0 67. 30 
2 5. 4176 0 18. O1 i. 2 5. 8490 0 &. 557 
$ 5. 4176 2. 258 196.8 5,3 5. 8490 2. 300 — 103. 0 
5. 4176 2. 258 196.8 5, 4 5. 8490 2. 300 — 103. 0 
l 5. 5003 0 257.1 4, 1 6. OO76 0 —134.4 
3 5. 5003 0. 5979 206, 0 4,3 6. 0O76 0.614 107. 00 
{ 5. 5003 5079 206, 0 4.4 6. 0O76 614 107. 00 
l 5. S732 0 262. 2 1,1 6.1474 0 137.0 
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wel. The lowest lving level belongs to T, and arises 
fom ),. Matrix elements for w-H7 that connect any 
igenvector of J), with any vector of the matrices of 
ye other determinants must vanish. Therefore 
ly levels with 7=1 are thus connected to T), and 
pence levels with 7=1, 3, and 4 belong to Ty. Sim- 
larly for all other levels. 

~The results in table I cover a wide range of situa- 
‘ons, and allow one to calculate magnetic properties 
quite readily. The method of presentation of the 
sults permits their adaptation in two ways. First, 
the eigenvalues can be adapted to any finite nonzero 
value of the spin-orbit coupling constant ¢. For posi- 
ve ¢ the problem corresponds to the octahedral case 
yith #2 configuration. (For negative ¢ the problem 
erresponds to the cubic case with /” configuration.) 
The adaptation is straightforward. Thus for D/¢ 





1i0-* and for the state N=2,3, W,/¢ 5.9820. 
W,/B 0.404, W.e/s? 12.4. Hence W, 
5.98206, W; 0.4048, W, 12.467/¢. If ¢=10° 
am, W, 5.9820 10° em™!, W, 0.4048, W, 
1248 107°. The second adaptation involves 


aterpolation of the results for values of D/¢ other 
than those covered here. This can be accomplished 
jy a graphical interpolation of the values in table I. 
These vary sufficiently slowly to permit interpolation 
f W., Wi, and W, for each state. 


3. Calculated Magnetic Susceptibility and 
Discussion of the Results 


The magnetic susceptibility per mole, Xv. of ions 
sith f? configuration was calculated [7] for a powder, 
a function of temperature and D). The results 
we plotted as x~' in figure 2. The results for our 
alculation for ¢=10° em™, and ¢=5 10? em” are 
riven for several values of J). Also some results of 
the theory of Penney and Schlapp [2] (effectivels 
©) are shown. For D/¢t lon!” 10-*, and 
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FIGURE 2. x7! versus temperature. 
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Gand Ds are taken from Penney and Schlapp theory (¢ , and crystal-field 


su tings comparable to C; and Dp, re spectively. 
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510°, with ¢= 10° em~', x"! is: practically inde- 
pendent of temperature up to 300°K, and equals 
5,102, 1,080, and 602, respectively. 

The matrix elements of H.+¢Z-S are given [2] in 
terms of the two undetermined parameters ¢ and D. 
The value of ¢ in a solid is uncertain, though the free 
ion value is often a useful first approximation. The 
correct value of ¢ for U’** lies probably [3] between 
600 em~' and 800 em. Thus the calculated sus- 
ceptibilities reported in this paper bracket the most 
likely value. Direct calculation, using eigenstates 
F?°H MsM,) obtained by standard methods [8] shows 
for the regular octahedron that the parameter 
dD y 2ee(r*)/(9ro°), indicating that D is negative. 

The calculations reported in this paper agree in a 
general way with the measurements and calculations 
of Hutchison and Candela [3] on [(CH,)4N], U Clg. 
They find a value of xy'=520, which is independent 
of temperature over the range measured, and explain 
this satisfactorily with their theory. Because the 
calculations reported in this paper always give a 
nondegenerate groundstate (T,), they also give a 
temperature independent X,;' for sufficiently large 
crystal-field splittings, and this for any value of ¢. 
The calculations reported in this paper shed some 
light on other measurements [9] of powder suscepti- 
bilities of compounds containing U**. These measure 
ments (on U(SOQ,). 4H,O, U(C.0,). 6H,O, UF,yNaF, 
(yJ)UFQNaF, and UFNaF) cannot be explained by 
a theory of the Penney and Schlapp type (i. e., 
(=o), nor by the present theory, though the present 
calculations somewhat improve agreement with the 
measurements on U(C,O,) 6H,O. This indicates 
that the finite value of ¢ does not play a very im- 
portant role in these systems, and that a more 
realistic crystal-field potential, involving terms of 
second, fourth, and sixth degree, is required for an 
understanding of these systems. 

The computation of the eigenvalues W and the 
eigenvectors was carried out on a 704 computer by 
Miss Irene Stegun of the Computation Section of the 
Applied Mathematics Division of the National 
Bureau of Standards. Her help was essential for 
this study, and is gratefully acknowledged. 
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